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FOREWORD 

The  work  described  in  this  report  was  performed  by  the  Virginia  Polytechnic 
Institute  and  State  University,  Engineering  Science  and  Mechanics  Department,  for  the 
Naval  Air  Development  Center  under  contract  number  N62269-85-C-0256.  The  principal 
investigator  was  Prof.  Norman  E.  Dowling.  DeRome  O,  Dunn,  Graduate  Research 
Assistant,  performed  important  portions  of  the  work.  Volume  II  was  taken  from  a  M  S 
thesis  written  by  K.  Ranganathan.  The  program  manager  for  NADC  was  L  W.  Gause. 
the  project  engineer  was  R.E.Vininc.  This  interim  report  covers 
work  that  was  perforr^ed  during  the  tiirie  period  October  1985  to  December  198 
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1.0  Introduction 


fatigue  crack  initiation  is  a  problem  which  causes  considerable  concern  since  it  occurs  in  man) 
materials  Most  fractures  in  machine  parts  arc  a  result  of  fatigue  cracking  I'atigue  failures  arc 
particularly  insidious  as  they  occur  at  nominal  or  average  stresses  which  arc  well  below  the  yield 
strength  of  the  matena]  and  hence  do  not  result  in  any  visible  signs  of  yielding  in  the  material 
f’ailurc  by  fatigue  begins  with  the  initiation  of  a  minute  crack  at  a  localized  spot,  usually  at  a  stress 
concentration  point,  and  gradually  propagates  through  the  section  of  an  engineering  component 
until  the  member  breaks 

One  aspect  of  material  behaviour  which  can  exert  a  significant  influence  on  the  fatigue  life,  but  has 
not  been  adequately  modelled,  is  the  relaxation  of  the  residual  stresses  present  in  the  material  Such 
relaxation  may  be  cycle  or  time  dependent.  A  number  of  methods,  both  destructive  and  non¬ 
destructive,  arc  available  for  measuring  residual  stresses  Hole  drilling  and  ring  coring  arc  methods 
that  arc  destructive  in  nature.  Some  non-destructive  techniques  used  for  measuring  residual  stresses 
include  acoustic  emission,  ultrasonic,  x-ray  diffraction  and  the  Barkhausen  noise  analysis  technique 
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Various  x-ray  diffraction  methods  have  been  used  in  the  past  for  different  applications.  However, 
measurements  base  yet  to  be  made  on  test  specimens  mounted  in  a  testing  machine  Measuring 
time  is  crucial  in  applications  which  require  the  study  of  time  dependent  relaxation  of  residual 
stresses 

Experiments  designed  to  study  the  time  dep>cndcnt  relaxation  of  residual  stresses  are  being  con¬ 
ducted  by  other  investigators  in  the  laboratory  on  specimens  containing  stress  raisers.  These  tests 
arc  being  done  on  notched  tensile  specimens  subjected  to  fatigue  loading 

The  objective  of  the  present  work  was  to  analylically  model  the  stale  of  stresses  dcscic'pcd  during 
fatigue  testing  of  notched  tensile  specimens  f’rom  this  model,  the  x-ray  response  to  the  strc'^s  slate 
of  the  specimen  was  determined  with  an  aim  of  measuring  the  x-ray  effeds  that  result  from  residual 
.stresses 


This  model  was  developed  in  two  parts.  In  the  first  part,  a  model  for  the  stress  distribution  on  the 
flat  or  plane  surface  of  a  notched  specimen  was  developed  for  the  region  between  the  notches  The 
corresponding  strain  distribution  was  then  calculated  from  Hooke's  law  of  elasticity.  These  calcu¬ 
lations  were  done  for  a  specimen  with  a  geometry  similar  to  that  used  in  experimental  work  This 
model  simulates  the  stress  state  which  develops  in  a  notched  tensile  specimen  during  experimental 
testing 

In  the  second  part,  the  x-ray  response  from  the  strain  distribution  was  first  determined  from  the 
variation  of  d,  the  lattice  spacing,  with  the  angle  il>.  'Ehis  variation  is  illustrated  as  plots  of  d  \  s 
sin^i//  Here,  i/<  is  the  angle  between  the  diffraction  vector  (which  is  the  bisector  of  the  incident  beam 
and  the  diffracted  beam)  and  the  normal  to  the  .sample  This  is  depicted  in  figure  1.  Ihe  second 
aspect  studied  in  this  model  was  the  relationship  between  the  stress  gradient  and  the  depth  of  pen¬ 
etration  of  the  x-ray  beam  Ihis  study  was  useful  in  estimating  the  error  involved  in  the  stress 
measurements  due  to  the  stress  gradient.  In  addition,  an  estimate  of  the  error  in  the  stress  values 
due  to  varying  radius  of  curvature  was  also  studied 
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An  important  part  of  this  simulation  was  the  determination  of  the  x-ray  effects  of  the  stress  state 
in  the  specimen  'ITiis  information  was  obtained  in  the  form  of  line  profiles  These  profiles  were 
simulated  for  the  individual  grains  on  the  flat  surface  of  the  specimen  in  the  region  between  the 
notches  using  a  finite  beam  size.  This  simulation  was  useful  in  the  study  of  the  effects  of  the  strain 
gradient  in  a  sp>ccimen  and  therefore  of  a  non-uniform  macrostrain.  The  theory  and  results  of  this 
simulation  arc  discussed  in  the  last  section. 
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2.0  Literature  Review 


Ilie  role  of  residual  stresses  in  the  premature  failure  of  metallic  components  has  been  recognized 
for  a  long  time.  There  is  a  great  need  for  the  mcaauremenl  of  these  stresses  for  a  number  o' reasons. 
These  include  stringent  quality  specifications,  the  reduction  of  load  safety  margins,  and  the  evolving 
application  of  fracture  mechanics  calculations.  Two  of  the  mo.st  prevalent  methods  |1|  used  >  r  the 
measurement  of  these  stresses  are  hole  drilling  and  ring  coring  These  two  methods  arc  sei.rde- 
structivc,  as  they  disturb  or  remove  material  from  the  component,  and  may  therefore  degrade  its 
seniccability.  The  methodologies  of  measurement  and  the  difficulties  associated  with  the  destruc¬ 
tive  techniques  necessitates  the  use  of  non  destructive  techniques  for  the  measurement  of  residual 
stresses. 

Numerous  non-destructive  techniques  have  been  used  for  measuring  stresses.  Some  of  these  meth¬ 
ods  are  -  x-ray  diffraction  |6),  ultrasonic  |2)  and  Darkhausen  noise  analysis  |3|.  Of  all  these  tech¬ 
niques,  x-ray  diffraction  is  the  only  non-destructive  technique  which  is  generally  reliable  and 
promising  for  a  wide  field  of  applications. 
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X-ray  diffraction  was  first  used  by  lister  and  Abom  |4|  in  (1Q25-1926).  I'he  diffracted  beams,  in 
those  days,  were  recorded  on  photographic  films  In  spite  of  the  great  advantage  of  being  nonde¬ 
structive,  x-ray  diffraction  has  not  been  used  w-idely,  until  recently.  This  was  due  to  the  excessive 
measurement  time  and  the  limitations  on  the  size  of  the  part  on  which  measurement  was  to  be  done 
(maximum  size  of  the  order  of  a  few  centimeters).  Today,  x-ray  diffraction  is  widely  used  due  to  the 
advancements  made  in  both  its  theoretical  and  technological  aspects 


2.  /  Residual  stresses 


Residual  stresses  are  defined  as  stresses  which  remain  in  a  material  after  all  applied  forces  arc  re¬ 
moved.  Considering  a  material  with  a  crystallographic  grain  structure,  such  as  metals  and  alloys, 
residual  stresses  may  be  classified  conventionally  into  three  types.  The  first  are  stresses  acting  over 
dimensions  of  the  order  of  a  few  grains.  These  stresses  arc  long  range  in  nature  and  may  be  termed 
macro-stresses.  Macro-stresses  may  develop  from  mechanical,  thermal  and  chemical  processes.  The 
second  kind  of  stresses  act  over  the  dimension  of  a  single  grain  and  are  termed  micro-stresses 
Micro-stresses  may  result  from  a  difference  in  the  mechanical  properties  of  the  different  phases  in 
a  material.  ITicy  are  manifest  as  localized  elastic  expansion  or  contraction  of  the  metal  lattice  may 
occur  as  a  consequence  of  processes  like  nitriding,  carburizing  or  phase  transformations.  The  third 
kind  of  stresses  range  over  a  few  inter-atomic  distances.  These  stresses  arise  from  the  varying  stress 
fields  of  individual  dislocations,  dislocation  pilcups,  kink  boundaries,  and  other  microstructural 
phenomena  which  are  of  a  discontinuous  nature.  The  macro  and  the  micro-stresses  are  measurable 
quantitatively  by  x-ray  diffraction  techniques.  They  may  be  detected  separately  by  x-ray  diffraction 
methods,  if  the  total  strains  from  all  phases  of  a  multiphase  solid  can  be  measured  The  volume 
fraction  of  all  the  phases  should  also  be  mea.surcd. 
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Surface  residual  stresses  exert  a  significant  influence  on  the  life  of  a  material  and  occur  at  notches 
and  scratches  Furthermore,  such  surface  defects  which  concentrate  additional  stresses  near  the 
surface  can  act  to  initiate  a  crack.  Compressive  residual  stresses  in  the  surface  are  beneficial  to 
components  prone  to  fatigue  failure  or  stress  corrosion  cracking,  as  these  stresses  have  to  be  over¬ 
come  by  the  applied  load  before  cracks  are  initiated.  On  the  other  hand,  residual  tensile  stresses 
are  usually  harmful  and  can  be  a  major  cause  of  brittle  fracture,  especially  in  fatigue.  Consequently, 
the  measurement  of  residual  stresses  is  of  considerable  importance  in  determining  the  load  carrying 
capabilities  of  the  material 


2.2  X-ray  diffraction 


The  principle  of  residual  stress  mca.suremcnts  by  x-ray  diffraction  is  based  on  the  fact  that  when  a 
polycrystalline  metal  or  ceramic  material  is  placed  under  stress,  the  elastic  strains  in  the  materia!  are 
manifest  as  a  change  in  the  lattice  spacings  of  its  grains.  The  stress,  be  it  externally  applied  or  res¬ 
idual  in  the  material,  is  taken  up  by  the  interatomic  spacings  X-ray  diffraction  techniques  can  ac¬ 
tually  measure  these  interatomic  spacings  which  are  indicative  of  the  macro  elastic  strain  in  the 
specimen.  The  corresponding  stress  values  are  obtained  from  these  elastic  strains  from  clastic 
stress-strain  relationships.  In  applying  these  relationships,  the  material  is  assumed  to  be  isotropic 
and  homogenous  in  all  of  the  volume  irradiated  by  the  x-rays. 

Monochromatic  x-ray  radiation  diffracts  according  to  Bragg's  law 

n).  —  2dun0  (2.0 1 ) 
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where*  ).  is  the  wavelength,  n  is  the  order  of  reflection,  and  d  is  the  interplanar  spacing  The 
diffraction  angle  20,  is  the  angle  between  the  diffracted  beam  and  the  incident  beam,  and  is  meas¬ 
ured  experimentally.  Figure  1  is  a  schematic  illustration  showing  how  x-ray  diffraction  occurs  and 
also  shows  the  effect  of  a  change  in  where  4>  is  the  angle  between  the  diffraction  vector  (which 
is  the  bisector  of  the  incident  beam  and  the  diffracted  beam)  and  the  normal  to  the  sample.  In  the 
figure,  a  and  /?  are  the  incident  and  diffracting  angles  respectively  If  planes  of  atoms  lying  parallel 
to  the  surface  |with  an  interplanar  spacing 'd')  satisfy  Bragg's  law,  diffraction  will  occur  as  shown 
in  Figure  1(a).  This  diffraction  occurs  from  a  thin  surface  layer  of  approximately  When  the 

specimen  is  tilled  with  respect  to  the  incident  beam  by  an  angle  i//,  a  different  set  of  grains  diffract, 
and  the  orientation  of  the  diffracting  planes  change.  W'hcn  a  uniform  stress  is  applied  on  the  sample, 
as  the  planes  that  are  diffracting  are  different,  the  value  of  d  is  changed  I  rom  Bragg  s  law,  a  change 
in  the  d-spacing  leads  to  a  corresponding  change  in  20.  ITic  presence  of  a  surface  .stress  in  the 
sample  is  reflected  as  a  shift  in  the  diffraction  peak  on  an  Intensity  -  20  plot  as  shown  in  Figure  1, 
Strain  in  the  sample  is  reflected  as  a  change  in  the  inter-planar  spacing  d,,,,  of  a  set  of  diffracting 
planes  (hkl).  This  variation  Ad  is  measured  from  Bragg's  law  as 

t  =  ■—  =  -  ( cot  0)  (2.02a) 

°hkl 


Solving  for  Af)  we  have, 


A0  =  — ( tan  0)  =  —  c  tan  ft  (2.02d) 

°hkl 

where  AO  is  the  shift  in  the  diffraction  peak.  From  the  value  of  the  20  peak  and  Bragg  s  law,  the 
spacing  of  the  planes  for  the  corresponding  orientation  can  be  obtained  From  equation  (2  02b). 
it  is  evident  that,  x-ray  diffraction  measurements  at  high  angles  have  a  greater  resolution  in  0.  This 
is  because  for  a  given  e,  as  0  increases,  tanO  increases  resulting  in  an  increased  sensitivity  and 
therefore  a  greater  resolution  in  0. 


T  he  theoretical  developments  necessary  for  the  treatment  of  lattice  strains  and  their  relation  to 
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stresses  follow  from  the  classical  theory  of  elastiniy.  The  state  of  stress  on  any  element  in  a  ho¬ 
mogeneous  isotropic  solid  (Figure  2a).  is  described  by  three  normal  components  o,,.  Oj,.  Oj,  and 
six  shear  components  <t,j.  0^,  oj, ,  o,,.  Oj,.  and  e,,.  The  stress  tensor  <t„  (i  =  1 .2,3;  j  =  1 ,2,3)  has  been 
proved  15|  to  be  a  second  rank  symmetric  tensor.  This  implies  that  <t„  =  a^,  ,  and  the  number  of 
independent  shear  stresses  is  reduced  to  three. 


The  relation  of  the  strain  equation  obtained  from  the  classical  theory  of  elasticity  to  x-ray  stress 
analysis  has  been  presented  by  Hillcy  16).  Figure  2b  illustrates  a  flat  metal  surface  wherein  0,1  , 
tTjj  ,  CT,,  and  t,|  ,  122  •  '■33  3re  the  principal  stresses  and  strains  respectively  For  any  stress  component 
<7,,,  from  elasticity  theory,  the  strain  in  any  direction  making  angles  <f)  and  t/'  with  <7,,  and  re¬ 
spectively  is  given  by 


[(<1, ,  cos^<^  +  '’12  +  '’22  ~  '’33I  sin^i/'  —  <7,^] 


-b  i-S*)  (^]  ]  -1-  ^22  "h  ^32^  T  -  2"  1 3  ^t)S  <j>  -f  ^23  5tn  <f>)  sin2\J'  (2.03) 

If  the  x-ray  diffraction  measurements  are  confined  to  the  surface,  then  plane  stress  conditions  arc 
assumed  to  exist  In  the  event  of  plane  stress  conditions  and  the  absence  of  shear  stresses,  equation 
(2.03)  reduces  to 


tS'2  2  2  2 

^H'  = - ^  '’22 I  ”22^ 


In  this  equation,  ^  is  the  spacing  of  the  planes  whose  normal  makes  angles  <t>  and  ij/  with  o,;  and 
<7,3  respectively.  S2/2  and  Si  arc  the  clastic  constants.  For  an  isotropic  medium, 

S2/2  =  (1  +  v)/F  and  S,  =  -v/F, 

where  F  is  the  Y'oung's  modulus  and  v  is  Poisson's  ratio.  When  \l/  =  ‘20"',  the  surface  stress  com¬ 


ponent  fT^  is  given  by 


2j.  ,  ■  2  , 

fff,  ~  O]  cos  -f-  (72  sin 
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Substituting  equation  (2.05)  into  equation  (2  0.5).  we  have, 

~ - 1 <'4.  S'n  ^  -  -^1  (oi  +  "2)  (2  06) 

“o  ^ 

When  ^  is  zero,  the  x-ray  beam  defects  the  strain  in  the  z  i.e,  the  c,,  direction.  This  strain  is  given 
by, 


*4,0  -  *^3.1  - 


(2.0^) 


Here  is  the  spacing  of  the  planes  reflecting  under  stress  at  ijf  =  0.  I  hc  strain  r,,  has  a  finite 
value,  given  by  the  Poisson  contractions  due  to  o,,  and  ojj  as 


®33  -  ~  •‘*1  ("11  +  ^22) 


(2.0,5) 


Subtraction  of  equation  (2.05)  from  equation  (2.06)  yields 


“  *33 


A 

2 


sin 


(2.09) 


This  is  the  basic  equation  which  relates  the  lattice  strains  to  stresses,  and  is  utilized  in  the  x-ray 
diffraction  method  of  stress  measurement.  Ihis  equation  is  for  conditions  of  plane  stress  I  x- 
pressing  the  strains  in  terms  of  plane  spacings.  we  have 


-d„ 


o 


A. 

E 


(I  v)  sin^^ 


(2.10) 


Equation  (2  10)  predicts  a  linear  variation  of  d  vs  sin^^  .  ITic  surface  stress  in  any  direction  on 
the  surface  of  the  specimen  can  be  obtained  from  the  slope  of  a  least-square  line  fitted  to  the  ex¬ 
perimental  data,  measured  at  different  As  multiple  -tilts  are  utilized,  this  procedure  is  known 
as  the  sinV  technique.  As  ^  „  and  d„  arc  very  nearly  equal  to  one  another,  (d^  ^  J  is  small 

compared  to  </„  .  The  unknown  spacing  </,  can  therefore  be  replaced  in  the  denominator  of  equation 
(2.10)  by  d^  ^  or  ,  .  For  most  materials,  it  has  been  reported  by  Noyan  and  C'ohen  |7|  that  the 
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error  introduced  by  the  elastic  strain  due  to  the  above  substitution  is  small  and  is  of  the  order  of 
0,1  ®/o.  Equation  (2.10)  can  now  be  written  as 


( 1  +  v)  sin  ^ 


(2.11) 


The  elastic  stress  constants  Sj/2  and  Si  are  known  to  vary  with  the  crystallographic  dircctitions  in 
an  anisotropic  material  As  the  strains  are  measured  in  particular  crystallographic  directions,  the 
mechanically  measured  values  of  these  constants  cannot  be  directly  applied  for  diffraction  meas¬ 
urements  It  is  hence  imperative  to  experimentally  determine  the  elastic  constants  in  the  measure¬ 
ment  of  stress  by  x-ray  diffraction.  The  Four-Point  Bend  Method  has  been  found  by  Prevey  (R| 
to  be  a  standard  and  a  simple  method  for  experimentally  determining  this  x-ray  elastic  constant 
parameter.  It  is  also  possible  to  calculate  them  from  single  crystal  clastic  constants  A  list  of  values 
of  the  elastic  constants  has  been  compiled  by  1  echnology  for  Energy  (Corporation  1^1  and  by  .lames 
listed  in  the  book  by  Noyan  and  Cohen  (7). 

The  manual  measurement  of  residual  stress  requires  the  accumulation  of  x-ray  counts,  for  different 
ij/  angles  with  the  first  ij/  angle  being  equal  to  zero.  The  accumulation  times  for  these  are  long.  If 
the  intensities  arc  collected  using  a  hardware  controlled  step  scan,  the  operator  time  arising  from 
long  collection  times  can  be  saved.  Thus,  with  a  system  that  includes  hardware  and  software  for 
on-line  computer  control  [lO],  the  operator  is  just  required  to  set  the  sample  in  the  diffractometer 
and  start  the  program  All  essential  procedures  like  peak  location,  selection  of  three  angles,  col¬ 
lection  of  x-ray  counts,  correction  of  intensities,  calculation  of  peak  angles  and  stress,  and  report 
of  the  data,  are  controlled  by  the  computer. 

The  last  decade  has  seen  an  increase  in  the  instrumentation  of  x-ray  diffraction,  offering  automation 
in  the  measurement  of  stresses  using  x-rays.  Instrumentation  such  as  Eastress  (1 1|  and  .sy.stems  sold 
by  Rigaku-Dcnki  of  Japan  offer  automation  nonetheless,  but  at  the  expense  of  precision  and  ver¬ 
satility,  Neither  method  measures  a  full  diffraction  peak.  Hence,  corrections  for  secondary  effects 
like  Ix)rcnt7,  polarisation,  sample  absorption,  beam  geometry  ,  resolution  of  the  K,  doublet  arising 
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out  of  anisotropic  cold  work,  etc.,  cannot  be  made.  7  he  computer  automated  stress  analysis  system 
reported  by  Koves  and  IIo  (1 2]  lack  on-line  control  for  locating  peaks,  which  requires  the  collection 
of  many  data  points  and  off-line  calculation  of  residual  stresses. 

The  use  of  position  sensiiive  detectors  makes  it  possible  to  record  an  entire  diffraction  peak  without 
moving  the  diffractometer  Position  sensitive  proportional  counters  (l.7|  not  only  process  the 
presence  of  a  photon,  but  also  its  location  along  a  line,  so  that  data  can  be  recorded  simultaneously 
from  a  range  of  ID  angles.  This  results  in  a  considerable  reduction  of  data  acquisition  time  Also, 
the  absence  of  the  motion  of  the  detector  to  record  a  peak  eliminates  the  heavy  and  expensive 
gearing  characteristics  of  a  conventional  diffractometer.  A  portable  x-ray  stress  analsvcr  based  on 
one  such  position  sensitive  detector  was  used  by  James  and  Cohen  (M]  for  residual  stress  meas¬ 
urements 

Various  methods  of  x-ray  diffraction  have  been  used  for  residual  stress  measurements  for  different 
applications.  However,  measurements  have  yet  to  be  made  on  test  specimens  mounted  in  a  testing 
machine.  Measurement  time  is  crucial  in  applications  which  require  the  .study  of  time  dependent 
relaxation  of  residual  stresses  Such  measurements  cannot  be  made  using  a  traditional  x-ray  system 
due  to  excessive  measurement  times  Moreover,  the  heavy  gearing  characteristics  of  such  a  system 
creates  a  problem  in  the  portability  of  the  equipment  to  the  test  site 

The  .stalc-of-thc-art  x-ray  stress  analysis  system  manufactured  by  Technology  for  Tnergy  Corpo¬ 
ration,  Knoxville,  TN,  is  being  used  for  making  such  measurements  on  notched  tensile  specimens 
This  system  is  based  on  the  work  of  James  and  Cohen  |15-IR]  The  system  consists  of  newly  de¬ 
signed  portable  diffractometer,  a  position-sensitive  proportional  counter,  signal -processing  elec¬ 
tronic  system,  an  x-ray  source  and  power  supply,  a  collimator  and  shutter,  an  articulated  arm 
diffractometer  support,  a  DFC  I^F-II/23  Plus  microcomputer  with  dual  floppy  disk  drives,  a 
printer-plotter,  a  graphics  display  terminal,  and  an  x-ray  safety  system  |9|.  All  components  arc 
mounted  in  an  air-conditioned  mobile  carl  An  entire  set  of  soflware  is  available  for  data  acquisi¬ 
tion.  calibration  and  data  analysis 


14 


NADC-88141-60  (Volume  II) 


A  distinct  advantage  of  the  advanced  x-ray  system  is  that  it  is  highly  compact  and  portable.  The 
instrumentation  can  be  mounted  for  work  on  a  specimen  in  a  mechanical  testing  machine  and  can 
remain  in  position  for  measurements  during  brief  pauses  in  a  fatigue  lest  ITiis  avoids  lost  data  due 
to  the  delay  that  would  otherwise  be  involved  in  removing  the  test  specimen  for  study.  It  also 
lessens  the  chance  of  accidental  damage  or  compromised  specimen  alignment,  both  of  which  arc 
more  likely  if  test  specimens  must  be  repeatedly  removed  and  remounted  in  the  testing  machine 
A  brief  description  of  the  experimental  set-up  used  for  making  measurements  along  the  edge  of  a 
notched  tensile  specimen,  at  the  bottom  of  the  notch,  is  given  by  Howling  and  Ranganathan  |19| 
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3.0  Stress  Distribution  Model 


3. 1  Introduction 


WTicn  a  polycnslallinc  metal  undergoes  deformation  with  the  strain  being  uniform  over  large  di";- 
tances  relative  to  the  grain  size,  the  elastic  portion  of  the  strain  causes  the  interplanar  spacings  to 
change  from  their  standard  stress  free  value.  ITic  change  is  dependent  on  the  magnitude  of  the 
applied  stress  This  change  is  indicated  by  a  shift  in  the  diffraction  lines  as  a  function  of  2P  position^ 
1  he  shift  in  the  Bragg  angle  position  (W)  due  to  the  strain  can  be  determined  experimentally  from 
which  the  clastic  strain  can  be  calculated.  Once  the  clastic  strain  is  known,  the  stress  can  be  calcu¬ 
lated  by  a  calibration  procedure,  or  from  stress-strain  relations 

In  the  analylical  work  done  in  the  present  work,  the  stress  state  in  a  notched  specimen  subjected 
to  uniaxial  loading  was  first  simulated,  which  is  presented  in  this  section 
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A  material  is  designated  as  elastic  if  it  obeys  Hooke's  law  of  elasticity.  In  this  limit,  the  ratio  of  the 
maximum  stress  to  the  nominal  stress  is  termed  as  the  stress  concentration  factor.  This  factor  is 
known  to  be  dependent  on  the  type  of  loading  and  the  geometry  of  the  specimen.  The  presence 
of  a  notch  leads  to  a  stress  variation  in  the  immediate  vicinity  of  the  notch.  At  large  distances  away 
from  the  notch,  the  stress  distribution  is  not  affected  by  the  notch;  for  example,  it  is  uniform  for 
simple  axial  loading.  The  shape  of  the  notch  in  a  notched  specimen  is  of  critical  importance  in  as¬ 
sessing  the  stress  gradients  present  in  the  vicinity  of  the  notch.  The  curvature  at  the  base  of  the 
notch  is  of  greater  importance  than  the  flank  angle.  In  the  case  of  a  double  edged  notch  specimen 
as  shown  in  Figure  .t,  the  ratio  of  half  the  distance  between  the  notches  'a",  and  the  radius  of  cur¬ 
vature  at  the  notch  tip  p.  is  a  non-dimensional  ratio  termed  as  the  "notch  cur\ature".  The  stress 
concentration  factor  has  been  found  to  depend  on  this  ratio  (2n|.  In  a  specimen  with  a  shallow 
notch,  the  stress  concentration  factor  is  independent  of  the  width  of  the  specimen  In  a  specimen 
with  a  deep  notch,  the  stress  concentration  factor  is  a  function  of  the  notch  curvature  alone,  the 
effect  of  the  depth  of  the  notch  being  insignificant  as  the  stresses  fall  off  rapidly  as  wc  go  farther 
away  from  the  notch.  As  a  limiting  condition,  when  h,  the  depth  of  the  notch,  is  zero,  the  stress 
(  oncentration  factor  becomes  equal  to  unity  . 

A  double  edged  notch  specimen  with  a  flank  angle  of  .10''  from  the  center  line  was  used  for  this 
investigation  1  he  notches  have  a  root  radius  of  p.  On  the  plane  surface,  the  notch  tips  are  separated 
by  a  distance  2a  A  schematic  illustration  of  the  specimen  geometry  is  shown  in  Figure  1.  Machine 
drawings  of  the  specimen  are  shown  in  Figure  4. 


3.2  Stress  distribution  calculations 


The  most  commonly  used  coordinate  .sy.sfcm  for  .stress  and  .strain  calculations  is  the  Cartesian  co¬ 
ordinate  system.  An  important  criterion  for  the  calculation  of  the  stresses  acting  on  the  surface  of 
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a  body  is  the  satisfaction  of  the  boundary  conditions.  In  most  cases,  as  the  body  is  curved,  it  is 
difficult  to  satisfy  the  conditions  using  the  Cartesian  system.  It  is  hence  preferable  to  use  a 
curvilinear  system  of  coordinates  for  such  calculations. 

The  three  orthogonal  directions  in  the  curvilinear  system  of  coordinates  are  u,  v  and  w.  These  and 
the  cartesian  coordinates  x,  y  and  z  are  assumed  to  be  differentiable  and  are  related  by  the  following 
relations: 


X  =  3' =  >>(h.v,w)  and  z  =  zfw.v.w)  (3.01) 

In  the  curs'ilinear  system  of  coordinatc.s,  the  surfaces  arc  defined  by  u  =  constant,  =  constant 
and  w  =  constant.  The  direction  of  u  at  any  point  is  given  by  the  normal  to  the  surface  u  - 
constant  Thus,  in  an  area  where  u  is  constant,  the  distance  between  two  adjacent  surfaces  sepa¬ 
rated  by  a  differential  du  is  given  by  A„du,  where  /i„  is  termed  as  the  factor  of  distortion  to  account 
for  the  curv'ilincar  distortion  resulting  from  differing  curvature  of  the  two  surfaces  The  factors  of 
distortion  for  the  surfaces  where  v  and  w  arc  constant,  are  pven  by,  h,  and  A,  respectively. 

The  direction  cosines  between  the  two  co-ordinate  systems  arc  given  by  the  cosine  of  the  angle 
between  the  linear  element  and  it's  corresponding  component  as  shown  in  Figure  5.  The  angles 
between  the  u  direction  in  the  curvilinear  system  and  the  x  direction  in  the  cartesian  system  is  given 
by 

cosIx-.m)  =  (3.02) 

du 

The  corresponding  angles  of  u  with  y  and  z  and  those  of  v  and  w  with  x,  y  and  z  can  be  calculated 
by  similar  equations. 

Considering  a  biaxial  system  such  as  a  bar  with  deep  hyperbolic  notches,  the  curvilinear  coordinates 
u  and  V  can  be  related  to  the  cartesian  coordinates  by  the  relations 
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jr  =  c  cosh  u  sin  v 
y  —  d  sink  u  cos  v 


(3.03) 


Thus  the  lines  u  =  constant  are  ellipses  and  the  lines  v  =  constant  are  hyperbolas;  c  and  d  are  nor¬ 
malizing  factors  between  the  two  coordinate  systems. 


Neuber  (20),  has  derived  the  stresses  in  cylindrical  coordinates  for  a  bar  with  infinitely  deep 
hyperbolic  notches,  in  uniaxial  tension. The  normal  stresses  in  the  u  and  v  directions  and  the  shear 
stress  arc  given  by  the  following  equations. 


A  2  2  2 

cosh  u  cos  v  (2/i  +  cos  —  cos  v) 
n 


A  2  2 

cosh  U  cos  V  (  cos  V  —  cos  Vp) 
h 


sink  u  sin  v  ( cos^v^  -  cos^v) 


(3.04) 


where  h  is  the  factor  of  distortion  and  is  found  from  the  Jacobian  of  the  two  systems  and  is  given 
by 

=  =  sinh^w  +  cos^',  (3.05) 


A  is  a  dimensionless  constant  given  in  terms  of  the  nominal  stress  p  (the  mean  stress  over  the 
narrowest  cross-section  under  tensile  loading)  as, 


A=p 


sin  >0 

-I-  sin  v„  cos  Vp  ’ 


(3.06) 


and  V,,  is  the  hyperbola  parameter.  A  definite  notch  contour  corrcsfxmds  to  every  value  of  v„  and 
this  is  related  to  the  notch  cunature  (-^)  by 


tan 


a 

r 


(3.07) 
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Using  the  two  relations  for  the  coordinate  transformation,  we  have,  for  a  point  on  the  hyperbola 
(and  hence  v  =  »„), 

=  (3,08) 

The  maximum  value  of  x,  for  each  value  of  y/d  is 

JCmax  =  c  cosh  u  sin  v 

i^)yld=0  =  ^  o 


Hence, 


(.I.OQ) 


In  order  to  express  the  parameter  v  as  a  function  of  x  and  y,  the  coordinate  transformation  relations 
are  used.  The  equation  of  the  hyperbola  is  pven  by. 


(x/c)^ 

.sin^v  cos^r 


(3,10) 


Solving  this  equation  for  v,  we  have. 


sin  V  = 


[  1  +  {xlcf  +  (yld)'^]  -  7[7+  (x/c)^  +  {yld)^Y  -  4{xlcf  1“''^ 

2 


(3.11) 


Note  that  the  positive  value  of  the  inner  square  root  yields  a  value  of  sin\v  greater  than  1  and  hence 
has  to  be  eliminated.  Thus,  summarizing  the  conditions  for  u  and  v,  we  have, 


u  =  sinh 


-I 


(y/0  ' 

COSV'o 


(3,12) 


V  =  Sin 


[  1  +  (xlcf  +  (yld)'^]  -  7[7+  ixlc)^  +  O'/rf)^]^  -  Mxirf 


-t/2 


(3.1.3) 
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The  values  of  the  stresses  in  curvilinear  coordinates  can  be  calculated  for  any  point  in  the  x-y  plane, 
in  the  region  between  the  notches. 

By  simple  coordinate  transformation,  the  corresponding  Cartesian  components  arc  found  to  be 


+  Oy  Oy  —  Oy 

Ox  =  ( - 2 - ^  ^ - 2 — 

-  <J„  0^—0^ 

: - )  -  ( - r - )rm  2n  +  It 


fT,  =  (  ■ 


(3.14) 


=  (  “  .  -  *  )fin  2-1  +  Tu^.eos  I't 


where  a  is  the  angle  between  the  tangent  to  the  curse  v  =  constant,  in  the  direction  u  increasing, 
and  the  x-axis.  and  is  given  by 


du 

sinh  u  sin  r 


(3.15) 


Hence, 


-1 

a  =  cos 


sink  u  sin  v 


] 


(3.16) 


Thus,  the  equations  presented  here  can  be  used  to  calculate  the  distribution  of  the  stresses  in  the  x 
and  y  direction  as  a  function  of  the  distance  from  the  notch  center  line. 
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3.3  Results  and  discussion 


The  stress  distribution  presented  thus  far  is  for  a  specimen  with  deep  hyperbolic  notches  in  an  in¬ 
finitely  wide  plate.  However,  in  order  to  extend  these  calculations  for  the  notch  specimen  with  a 
geometry  as  shown  in  Figure  3  ,  the  following  approximations  have  to  be  valid, 

•  F.xtension  of  the  calculations  for  a  specimen  with  infinite  geometry  to  a  specimen  with  finite 
geometry. 

•  Approximation  of  a  noteh  with  (lank  angles,  to  a  hyperbolie  notch 
Justifications  for  the  above  two  approximations  follow. 

The  stress  concentration  factors  for  two  different  notch  curvatures  arc  shown  in  1  able  1 .  The  values 
listed  are  for  the  cases  of  a  specimen  with  finite  and  infinite  width  respectively.  These  values  have 
been  established  from  the  Handbook  of  Stress  Concentration  Factors  pi).  The  stress  concentration 
factor  K,  was  first  obtained  for  a  flat  tension  bar  with  opposite  U-shaped  notches  for  the  a/^i  and 
w/a  which  is  pertinent  to  our  specimen  geometry,  where  w  is  the  width  of  the  specimen  The  stress 
concentration  factor  for  the  corresponding  V-shaped  notch  with  the  (lank  angle  a,  was  then  es¬ 
timated  from  a  plot  of  K„  vs  K,  for  different  flank  angles  o,  obtained  from  the  same  handbook. 
As  is  evident  from  the  tabulated  values,  the  error  involved  in  extending  the  calculations  from  an 
infinite  model  to  one  of  finite  width  is  approximately  3  percent. 

Kikukawa  |22I  conducted  tests  on  a  photoclastic  model  made  with  H-shaped  notches  having  a'p 
=  3.0  and  hjp  ranging  from  10  0  to  0.25.  He  reported  that  the  largest  increase  in  the  stress  con¬ 
centration  factor  for  notches  of  finite  depth  as  compared  to  infinitely  deep  notches  was  7  percent. 
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Tabic  I 


Stress  concentration  factors  for  finite  and  an  infinite  model. 


K 

tip  =  2.0 

tip  -  1.0 

Infinite 

Model 

1. 96 

1.5.‘i 

Finite 

Model 

2.02 

1.61 

Percentage 

Firror 

.10 
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Figure  6.  Neuber's  estimate  of  stress-concentration  factors  for  a 
tensile  bar  Kith  U-shaped  notches  (Fbnn  Roll  [24]). 
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Neuber  |201  made  an  analyliral  estimate  of  the  stress  concentration  factor  for  a  specimen  with 
symmetric  li-shapcd  notches  having  a'p  =  3.0  and  h'p  ranging  from  0  to  4,0.  Figure  6  shows  this 
estimate  of  the  stress  concentration  factor.  From  the  figure,  it  can  be  seen  that  the  change  in  the 
stress  concentration  factor  between  h/p  =  10  and  2.0  is  approximately  8  percent.  A  point  worth 
noting  in  this  regard  is  that  the  stress  concentration  factor  obtained  by  Kikukawa  (23|  for  infinitely 
deep  U-shaped  notches  was  the  same  as  the  value  obtained  for  infinitely  deep  hyperbolic  notchc*; 
In  all  cases,  Kikukawa  |231,  as  reported  by  Flynn  and  Roll  |24|,  found  that  for  certain  a/p  ratios, 
the  maximum  stress  concentration  factor  occurs  with  notches  of  finite  depth  rather  than  with  infi¬ 
nitely  deep  notches  Tbe  difTcrcnce  in  the  stress  concentration  factors  for  the  two  cases  is  not  sig¬ 
nificant  The  stress  distribution  in  the  region  between  the  notches  is  dependent  on  the  stress 
concentration  factor,  with  the  maximum  values  of  the  stress  being  in  the  vicinity  of  the  notches 
As  the  difference  in  the  stress  concentration  factors  for  the  two  models  is  not  significant,  the  stress 
distributions  for  the  two  models  are  similar.  Hence,  extension  of  the  calculation  for  a  specimen  with 
infinite  geometry  to  a  specimen  with  finite  geometry  is  a  valid  approximation 

The  stress  distribution  obtained  for  o,  and  o,  from  equation  (3.14)  the  calculations  are  plotted  as 
a  function  of  the  distance  from  the  centerline  in  Figure  7.  This  stress  variation  is  for  the  surface 
A  of  the  specimen.  Figure  7  illustrates  the  variation  of  both  the  x  and  y  component  of  the  stress 
ie.,  (T,  and  in  the  region  between  the  notches  on  surface  A  of  the  specimen  The  co-ordinate 
axis  is  as  shown  in  Figure  3,  with  the  ongin  nemg  at  the  centerline  on  surface  A  ,  in  the  region 
between  the  notches.  The  stresses  were  plotted  for  various  heights  (in  inches),  above  the  x-axis  ie., 
for  four  different  values  of  y. 

The  variation  in  e>^  is  of  importance  from  the  point  of  view  of  the  slit  geometry  being  used  in  the 
experimental  x-ray  measurements.  A  0.5mm  by  5mm  (0.02'  by  0.2")  rectangular  slit  is  one  of  the 
slit  geometries  used  in  the  experiments  being  conducted  at  Virginia  Tech  on  the  I  F.C  x-ray  system 
bv  other  in'  cstigators.  A  cross-section  of  this  slit  is  shown  in  Figure  7.  The  5mm  dimension  in  the 
longitudinal  direction  corresponds  to  0.1  inches  above  the  centerline.  Fhc  variation  in  the  peak 
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Distance  from  the  centerline  (inches) 

itrcss  distribution  in  t!>c  region  between  tlie  notches  for  diflerent  lieights  above 
le  centerline. 
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stres'.  corre<iponding  to  y  =  0.0  and  >  =0.1  is  approximately  5''o  At  regions  farther  away  from  the 
notch,  this  variation  is  seen  to  he  negligible.  This  implies  that  the  above  orientation  of  the  slit  is 
good  for  the  x-ray  measurements.  Also,  the  negligible  variation  in  the  stress  up  to  a  certain  height 
above  the  x-axis  helps  in  justifying  the  approximation  of  a  V-notch  with  a  flank  angle  to  a 
hyperbolic  notch  This  approximation  is  however  not  valid  beyond  a  certain  height  above  the  notch 
centerline,  due  to  the  considerable  change  in  the  peak  stress  as  seen  in  the  dotted  curves  of  Figure 
7.  Any  misalignment  in  the  slit  orientation  along  the  y  axis  would  cause  detection  of  the  variation 
in  the  peak  stress  corresponding  to  y  greater  than  0  1 . 

The  stress  distribution  obtained  so  far  is  used  as  the  basis  for  the  x  r.ay  diffraction  model  dealt  with 
in  the  next  section 


3.4  Calculation  of  strains 


As  strain  is  the  quantity  which  is  measured  experimentally,  it  is  important  to  establish  relations  for 
strains  in  the  three  orthogonal  directions.  This  is  done  using  Hooke's  law  of  clasticitv  '1  he  cor¬ 
responding  elastic  strains  are  obtained  from  the  following  equations  In  a  plane  stress  condition 

fT,  =  0. 

rv  =  y(n,-vn,) 

C;  =  - n,.) 

where  v  and  F  are  the  Poisson's  ratio  and  the  Young's  modulus  respectnely  V  alues  for  ^  and  1 
for  74’’.^  A1  were  obtained  from  the  book  by  Byars  and  Snyiler  |25|  and  are  given  in  I  able  2  Note 
that,  although  rr,  =  0,  there  is  a  strain  in  the  7  direction  (r,)  due  to  Poisson  effect 
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distribution  in  the  region  between  the  notches  as  a  function  of  distance 
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The  plot  of  the  strains  as  a  function  of  the  distance  from  the  centerline  is  shown  in  Figure  8  The 
strains  in  the  x  and  z  directions  are  compressive  due  to  Poisson  contraction.  The  strains  thus  ob¬ 
tained  were  used  in  determining  the  d-spacing  profile  from  which  the  intensity  profiles  were  simu¬ 
lated 
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4.0  X-ray  Diffraction  Model 


The  d-spacing  profiles,  calculated  from  the  strain  distribution  described  in  the  previous  chapter,  arc 
presented  in  this  section.  An  analytical  model  which  simulates  the  state  of  stress  developed  during 
experiments  conducted  on  notched  tensile  specimens  is  described.  The  x-ray  response  from  this 
model  was  obtained  in  order  to  control,  study  and  interpret  the  x-ray  effects  that  arc  seen  during 
testing.  The  presence  of  a  stress  in  the  specimen  is  indicated  by  changes  in  the  d-spacings  of  the 
diffracting  planes  The  line  profiles  resulting  from  the  d-spacing  variation  on  the  flat  surface  (A)  of 
the  specimens  were  simulated.  These  profiles  were  obtained,  as  x-ray  effects  from  individual  grains, 
for  a  beam  of  finite  si^e.  They  arise  as  a  re.sult  of  a  strain  gradient  present  at  the  specimen  surface 
The  effects  of  the  depth  of  penetration  of  the  x-ray  beam  on  the  stress  gradient  occuring  on  surface 
B  (edge)  of  the  specimen  is  also  discussed  in  this  section.  Sample  calculations  for  a  7475  Al 
notched  tensile  specimen  arc  also  presented. 


35 


NADC-881 41-60  (Volume  II) 


4. 1  Model  for  residual  stress  measurement 


The  strains  obtained  from  the  stress  distribution  model  were  used  to  calculate  the  variation  in  the 
inter-planar  spacings.  From  the  x-ray  defmition  of  strain,  it  is  known  that  the  shift  of  the  diffraction 
lines  to  new  20  positions,  i.e.,  the  difference  in  the  d-spacings  for  the  stressed  and  the  unstressed 
state  is  related  to  the  uniform  macro-strain  that  causes  the  shift  by  the  relation. 


il'  - 


(4.01) 


where, 

*  is  the  strain  in  a  direction  making  angles  <f>  and  with  n,,  and  n,,  directions  (Figure  2), 

^  is  the  spacing  under  stress  of  the  planes  normal  to  the  direction  making  angles  and  4'  "ith 
and  ffj,  respectively,  and 

is  the  spacing  of  the  same  planes  in  the  absence  of  stress. 


When  =  0,  the  x-ray  beam  detects  the  strain  in  the  z  direction.  I  his  strain  is  given  by 


(4.02) 


Here  ,  is  the  spacing  of  the  planes  reflecting  under  stress  at  ^  =  0.  The  relation  for  the  difference 
between  the  two  strains  is  given  by  equation  (2.09)  as. 


*’4'  2 


(4.03) 


where  is  the  surface  stress  component.  Again,  from  equation  (2.10)  we  have, 


— - — . - - - —  + 


(4.04) 
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The  above  equation  is  based  on  the  following  assumptions. 

•  The  material  is  isotropic  and  homogeneous  in  all  of  the  volume  irradiated  by  x-rays. 

•  Since  measurements  are  confmed  to  a  small  depth  below  the  surface,  plane  stress  conditions 
are  assumed  to  exist. 

•  Shear  stresses  are  absent  in  the  irradiated  volume. 

The  stress  can  be  obtained  from  the  slope  of  a  least-squares  line  fitted  to  experimental  data 
measured  at  various  \}/  if  the  constants  F,,  and  »  are  known  The  lattice  spacing  of  the  planes  in 
the  unstressed  state,  ,  may  not  be  available  in  practice.  1  he  lattice  spacing  measured  at  i/-  =  0 
is  substituted  for  It  has  been  reported  by  Noyan  and  Cohen  (7|  that  for  an  elastic  strain,  the  error 
introduced  due  to  this  substitution  is  only  0. 1  percent.  Fquation  (4.04)  can  now  be  written  as 

- — j— —  (4  05) 

(1  -1-  v)  sin  \li  ^4.0 

For  the  purpose  of  this  simulation,  the  value  of  rf,  was  used  as  a  .scaling  factor  in  the  calculation 
of  the  d-spacing  profile.  The  values  of  d,  for  reflections  from  the  (200),  (400)  and  (3.1.1)  planes  were 
obtained  from  Pciirson's  Handbook  (26j  These  values  and  the  wavelength  of  the  radiation  used  arc 
listed  in  Table  2  Table  2  also  contains  the  values  of  the  Young's  modulus  F  and  Poisson's  ratio 
V  used  in  the  calculations.  All  the  parameters  arc  listed  for  a  7475  A1  notched  tensile  specimen 

The  fundamental  equations  in  the  x-ray  diffraction  measurement  of  stresses  arc  now  used  to  develop 
a  model  for  the  determination  of  residual  stresses  on  the  flat  face,  surface  A  (Figure  3)  of  the  spec¬ 
imen  in  the  region  between  the  notches  and,  on  the  edge,  surface  B  of  the  specimen  at  the  bottom 
of  the  notch  (notch  (ip).  This  model  will  be  of  help  in  the  analysis  of  the  .stress  .state  rc.sulting  from 
the  experiments  being  conducted  in  the  notched  specimen.  fYom  equation  (2  03),  the  strain  in  any 
direction  can  be  determined.  A  coordinate  system  in  which,  the  unprimed  quantities  refer  to  the 
specimen  axi";  and  the  primed  quantities  imply  laboratory  coordinates  is  used,  following  the  con- 
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Table  2.  Material  parameters. 


Item 

7475  Al 

Radiation 

Cu 

(A) 

1.54178 

Rectangular  * 

.5  by  5.0 

Slit  Size,  mm 

do  (A) 

2.0235 

(400) 

l.Oll? 

(300) 

0.7788 

F,  (ksi) 

1 0300 

V 

0..33 

*  5mm  dimension  parallel  to  specimen  thickness. 
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vcntion  established  by  Dolle  |271.  The  strain  e,,.  in  the  direction  normal  to  the  diffracting  planes  can 
be  determined  from 


(4.06) 


d-  (til  cos^((>  +  tij  sin^^"  +  *22  sin^<^  —  C33)  sin^i(f  +  (C13  cos <(>  +  033  sin  sin  2</'] 

Assuming  the  ab.sence  of  shear  strains,  we  have 

*12  =  *1.1  =  C23  =  0 

Using  this  condition,  equation  (4.06)  becomes 


£33, (<t, »/-)=- 

=^£33  +  (til  co5^(()  4*  £32  4>  —  C33)  sin  i// J 


(4.07) 


and  hence  can  be  expressed  in  terms  of  the  strains  by, 

+  '33') 

+  *33)  +  (*n  cos^ <t>  +  £32  sin^d>  ~  *33)  sin^ii'] 


(4.0R) 


Note  that  £'3,(<^,  0)  =  £33 

The  d-spacings  at  each  point  along  the  x-axis  can  be  calculated  from  equation  (4.0R)  Fqualion 
(4.08)  relates  to  *  to  principal  strains.  The  strains  calculated  in  the  last  chapter  (Figure  8)  were 
in  the  x.  y,  7.  system.  The  principal  directions  for  cases  A  and  B  arc  different.  In  order  to  use 
equation  (4.08)  for  both  cases,  two  orthogonal  coordinate  systems  were  chosen  in  a  principal  strain 
space  i  c  ,  c,,.  c^j  and  £,,  space  which  were  then  related  to  the  .strains  c,  ,  £,  and  t,  The  values  of 
the  strains  were  obtained  from  the  strain  distribution  in  Figure  8.  A  figure  of  the  specimen  with  the 
corresponding  coordinate  systems  for  surfaces  A  and  B,  is  shown  in  f  igure  9. 
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Figure  9. 


(b) 


(a) 


A  figure  of  the  specimen  showing  the  coordinate  systems. 
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4.1.1  Plane  surface  (surface  A) 

With  reference  to  the  co-ordinate  system  in  Figure  9(a),  we  have  for  measurements  on  surface  A, 

®ii  = 

'22  =  Cy 

'32  =  h 


Substituting  these  into  equation  (4.07)  we  have. 


'33 


=[c^  -t-  -f  tj.  sin^d>  -  '.?)  sin^i//  ] 


(4.09) 


It  is  easy  to  verify  that, 


<^=  0” 

,  ^  =  0' 

•  '33'  =  '7 

</.  =  0” 

,  f  =  90" 

■  '33-  =  'x 

(4.10) 

<t>  =  90' 

,  (^  =  90" 

,  Cjy  =  ty 

as  is  required.  The  three  different  sets  of  values  of  and  correspond  to  three  different  positions 
of  the  x-ray  source  with  respect  to  the  specimen  surface.  These  three  positions  correspond  to  the 
measurement  of  the  d-spacings  in  each  of  the  three  orthogonal  directions 


The  d-spacing  in  any  direction  4>.  ij/  an  surface  A  is  obtained  from  equation  (4  09)  as, 

A 

+  '7)  +  ('s  cos^(/.  +  ty  sin^d-  -  tp  sin^i//]  (4.1 1) 

This  equation  was  used  to  calculate  the  d-spacing  for  ,^  =  0"  and  90"  at  seven  i/'  values.  These 
values  arc  tabulated  in  columns  II  and  IV  of  lablc  .1, 
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Table  y.  d-spacing  as  a  function  of  tl>  at  (-0.25,  0,  0) 


sin^  ^ 

<* 

0.00 

0.77535 

0.77535 

0.77535 

0.77535 

0.07 

0.77535 

0.77628 

0.77628 

0.77535 

0.25 

0.775.35 

0.77885 

0.77885 

0.77535 

o..5n 

0.77535 

0.782.14 

0.782,34 

0.775.35 

0.75 

0.77535 

0.78.584 

0.78.584 

0.775.35 

0.9,1 

0.77535 

0.78840 

0.78840 

0.775.35 

1  00 

0.77535 

0.789.34 

0.789.34 

0.77535 
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4.1.2  Edge  (surface  B) 

From  the  co-ordinate  system  shown  in  Figure  9(b),  for  measurements  on  surface  B, 

c,,  = 

Cjj  =  r~ 

*^3.1  = 


Consequently,  equation  (4.07)  now  becomes. 


~[*x  +  (''-y  4>  4-  sin^<^  —  tji)  sin^</'] 


(4.12) 


It  can  again  be  verified  for  consistency  that, 


,^=  O' 

.  ^  =  0”  , 

'  ^3r  ~ 

<^  =  0' 

o 

11 

,  tjy  =  ’^y 

(4.13) 

(t>  =  90” 

.  i  =  90" 

■  '^33'  ” 

Thus,  the  d-spacing  in  any  direction  4>,  ^  on  surface  B  can  be  obtained  as, 

p 

+  'x)  +  (c^  cos^if.  +  c,sin^</>  -  C;j)  sin^i/']  (4.14) 

This  equation  was  used  *o  calculate  the  d-spacing  for  <f>  ~  0"  and  90'  at  .seven  ^  values  These 
values  are  tabulated  in  columns  III  and  V  of  Tabic  .7. 
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4.2  Discussion  of  the  d  r5  sin^  ij/  results 


The  d-spacing  values  calculated  from  equations  (4.11)  and  (4.14)  are  tabulated  in  Table  ?>. 
Equations  (4.1 1)  and  (4.14)  are  dependent  on  the  strains  in  the  three  orthogonal  directions  x,  y  and 
z.  As  the  strain  distribution  (Figure  8)  is  a  function  of  position  along  the  x-axis,  the  d-spacing  values 
can  be  calculated  from  equations  (4.11)  and  (4.14)  for  various  positions  along  the  x-axis.  The  values 
listed  in  the  table  were  calculated  from  the  strain  value  at  x  =  -0.25  (Figure  8)  This  corresponds 
to  a  point  at  the  bottom  of  the  notch. 

Equations  (4.1 1)  and  (4.14)  can  be  utilized  to  obtain  the  d-spacing  (x.v )  at  any  pc'int  on  the 
surface  of  the  sample.  As  diffraction  takes  place  from  a  thin  surface  layer  (~  20//m)  plane  stress 
conditions  were  assumed  in  the  above  calculations.  The  geometry  of  the  specimen  chosen  for  this 
simulation  (Figure  3)  also  justifies  the  assumption  of  plane  stress,  as  the  notch  radius  is  not  small 
compared  to  the  thickness  of  the  specimen.  F^ch  set  of  values  of  <l>  and  4'  correspond  to  meas¬ 
urement  of  the  d-spacings  in  each  of  the  three  orthogonal  directions.  Consider  for  example,  the 
d-spacing  value  for  =  0"  and  4/  =  90'  on  surface  A.  This  value  is  a  result  of  the  strain  in  the  x 
direction  From  equation  (4  13),  it  can  be  inferred  that  the  d-spacing  measured  on  surface  B  at  (fi 
=  0'  and  4‘  -  O'  is  a  result  of  the  strain  t,.  As  the  same  point  and  thus  the  same  set  of  planes  arc 
being  measured  on  both  surfaces  A  and  B,  the  d-spacing  resulting  from  the  strain  (r.,  in  this  case) 
ought  to  be  the  same  This  is  found  to  be  true  by  comparing  the  values  in  Table  3  This  then 
validates  the  above  calculations  The  results  can  be  summarized  as  follows. 


/I  A 

4),o  =  n 

fl  fl 

4n.n  =  ^0,0 


(4.15) 


where  the  superscripts  'A'  and  B'  correspond  to  surfaces  A  and  B  respectively.  This  implies  that 
when  4'  ~  ^  '''  independent  of  <i>  which  is  seen  to  be  true  from  equation"!  (4.13)  and  (4  14).  Also. 
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A  B 

4), 90  =  4),n 

A  B 

^0,90  =  90  (4  If') 

A  B 

4d,0  =  *^0,90 

This  equation  implies  that  the  same  set  of  planes  arc  being  measured  in  both  cases  for  the  limiting 
conditions.  The  results  of  Table  3  are  plotted  in  Figure  10  Fhc  different  nature  of  the  strains  such 
as,  ty  being  tensile  and  c,  being  compressive  leads  to  varying  slopes  of  the  two  lines  in  f  igurc  10. 

A  f 

From  Figure  10  it  can  be  seen  that  the  d  vs  sin^i//  plots  arc  colinear  for  and  .  Similarly,  thc> 

f  A 

are  colinear  for  4).#  and  4o,^  This  is  due  to  the  fact  that  r,  and  r.,  arc  equal  at  the  bottom  of  the 
notch  i.e.,  at  x  =  -0.25  in  Figure  8.  It  can  therefore  be  inferred  that  the  contribution  of  a  particular 
strain  component  plays  an  important  role  in  the  determination  of  the  intcrplanar  spacings.  Ibis  is 
seen  to  be  true  from  equations  (4.11)  and  (4.14).  As  the  sensitivity  of  the  experiment  is  ~  0  00002 
A.  the  effects  illustrated  in  Figure  10  can  be  detected.  In  this  context  it  is  worth  noting  that  attaining 
precision  in  the  x-y  positioning  of  the  x-ray  beam  on  the  specimen,  particularly  at  the  bottom  of 
the  notch,  is  of  prime  importance. 


Depth  of  penetration  calculations 


X-ray  diffraction  techniques  are  confined  to  surface  measurements  Ibc  strong  attenuation  due  to 
absorption  limits  the  depth  of  penetration  of  fhc  x-ray  beam,  so  that  diffraction  occurs  from  a  thin 
surface  layer  of  the  material  The  presence  of  a  stress  gradient  in  the  near  surface  layers  of  the 
sample  lead  to  a  non-linearity  in  the  d  vs  sin’^  plots  discussed  in  the  last  section.  Tbis  section  deals 
with  the  effect  of  the  depth  of  penetration  on  the  gradient  occuring  on  .surface  B  of  the  specimen 
i.e,  for  edge  measurements.  To  begin  with,  a  review  of  the  theory  of  absorption  is  given 
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4.4  Theory  of  absorption 


In  a  reflection  technique  by  x-ray  diffraction,  the  irradiated  volume  that  is  observed  includes  the 
attenuation  effect  due  to  the  absorption  of  the  x-rays  in  the  material.  Consider  an  x-ray  beam  of 
cross-sectional  area  Ao,  irradiating  the  surface  of  a  sample  at  an  angle  of  incidence  0.  The  irradiated 
surface  area  A,  is  then  given  by, 

(4.17) 


Figure  1 1  illustrates  such  a  region.  The  maximum  absorption  is  dependent  on  the  linear  absorption 
coefficient  //  of  the  material.  A  number  of  sampling  planes  of  area  A,  at  various  distances  can  be 
imagined.  Sample  absorption  reduces  the  intensify  by  a  factor 


exp 


where,  ^  is  the  linear  absorption  coefficient  of  the  material  and  Z  is  the  depth  at  which  the  ab 
sorption  is  considered.  The  differential  intensity  is  then  given  by 


-2nz 
sin  0 


d? 


(4.18) 


where,  /,  is  the  incident  intensity. 

If  a  semi-infinite  slab  is  considered,  the  total  diffracted  intensity  is  obtained  by  integrating  equation 
(4.18)  over  an  infinite  thickness.  This  intensity  can  be  found  to  be 

/  =  -^  (4.19) 

2/1 
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Here,  /,,  and  ^  are  constant  for  all  reflections  and  is  independent  of  0.  ITic  absorption  factor, 
1, 2/(,  for  a  scmi-infiiiilc  slab  making  equal  angles  with  the  incident  and  the  diflractcd  beam,  is  hence 
found  to  be  independent  of  6.  This  is  due  to  the  balancing  of  two  opposing  effects.  At  small  0. 
while  the  specimen  area  irradiated  by  an  incident  beam  of  fixed  cross-section  is  large,  the  effective 
depth  of  penetration  is  small.  At  large  0,  the  opposite  is  true.  The  above  equations  were  obtained 
from  CuUity  [2')]. 


When  the  specimen  is  tilted  by  an  angle  the  resulting  absorption  factor  is  given  |2S|  by  (1  - 
tarn/'  cotf^).  The  corresponding  equation  for  intensity  is 


r//^  = 


-f  Ip) 


.  expj^  -  /ir[  - 


1 


1 


sin(0  +  iji)  sin(fl  —  il>) 


dr 


(4.201 


Effects  due  to  this  absorption  and  factors  like  lT>rcnl7,  polarization  have  been  taken  care  of  in  the 
software  for  the  T.E.C.  instrument  (9|. 


The  total  intensity  diffracted  by  a  slab  of  material  between  this  layer  and  the  surface,  expressed  as 
a  fraction  of  the  total  diffracted  intensity  is  pven  I?)  by 


—  exp 


[1,1]' 

-  /‘Zl 

sin(0  +  f )  sin(fl  -  ip)  \ 

(4  21) 


The  effective  depth  of  penetration  maybe  defined  as  the  thickness  that  contributes  99%  of  the 
diffracted  intensity.  It  is  related  to  the  linear  absorption  coefficient  ;j,  ip  and  the  Bragg  angle  0. 
Tbus,  from  the  above  equation,  for  a  given  0.  /i  and  \p,  when  (7,  =  0.99,  z,  =  where  S  is  the 
depth  of  penetration. 


I  hc  effect  of  the  gradient  on  the  depth  of  penetration  on  the  stress  gradient  is  studied  for  meas¬ 
urements  on  the  edge  (surface  B)  of  the  specimen  As  there  is  not  expected  to  he  a  gradient  in  the 


49 


NADC-881 41-60  (Volume  II) 


stress  normal  to  surface  A,  the  effect  of  the  depth  of  penetration  for  measurements  on  surface  A  is 
not  important. 


4.4.1  Effect  of  the  depth  of  penetration  of  x-rays 


The  arrangement  of  the  specimen  for  edge  measurement  at  the  bottom  of  the  notch  is  shown  in 
Figure  12.  The  values  of  the  depth  of  penetration  calculated  for  ^  =  O',  where  it  is  greatest,  arc 
listed  in  T  abic  4.  This  table  also  contains  the  values  of  the  linear  absorption  cocfftcient  in  aluminum 
and  titanium  for  two  x-ray  radiations  1291-  The  value  of  0  used  in  equation  (4.21)  for  the  calculation 
of  S  is  also  given  in  the  table. 

With  reference  to  Figure  7,  and  considering  the  greatest  depth  of  penetration  listed  in  Table  4,  the 
decrease  in  the  stress  within  the  penetration  depth  is  not  large.  Due  to  this  stress  gradient  in  the 
surface  layers  of  the  sample,  the  diffraction  angles  that  are  measured  represent  samplings  of  different 
stress  values,  affecting  the  residual  stress  results  by  a  small  amount.  For  an  infinitesimal  x-ray  beam 
incident  at  the  bottom  of  the  notch,  as  in  Fig.  12,  and  for  the  slit  geometry  being  used  on  the  TF,(' 
x-ray  system,  the  errors  encountered  in  the  stress  measurement  by  x-rays  due  to  the  stress  gradient 
effect  do  not  exceed  about  2%  of  the  measured  stress  value. 

Another  factor  which  affects  the  accuracy  in  stress  measurements  using  x-ray  techniques  is  the  ra¬ 
dius  of  curs  aturc  of  the  sample.  This  effect  was  studied  for  two  different  notch  radii,  for  both  alu¬ 
minum  and  titanium  specimens.  The  results  of  this  study  arc  shown  in  Appendix  A.l.  The  error  in 
the  stress  values  due  to  varying  radius  of  curvature  was  found  to  be  of  the  order  of  2  ksi.  for  any 
value  of  measured  stress.  For  the  counting  statistics  available  on  the  1  Ii(,'  system,  this  error  is  of 
a  very  small  magnitude.  Hence,  this  effect  is  not  of  importance  for  measurements  made  using  the 
TFC  x-ray  system 
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Figure  12 .  Arrangement  for  measurements  on  the  edge  of  the  specimen 
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Table  4.  ^  alues  of  the  depth  of  penetration  in  Al  and  I'i  for  two  dilTcrent  radiations  used. 


(a)  Radiation  :  CuK,  at  V'  =  0° 


Materials 

Planes 

20,  deg 

fx,  cm  ' 

6,  10  ’ in, 

for  absorption  of; 

95% 

99  "i, 

Al 

f511/.W) 

162.5 

131.5 

4.4.3 

6R1 

a  —  hep 

1i 

(213) 

139.5 

926.2 

0.597 

0  91R 

(b)  Radiation:  CrK.  at  i/i  =  0° 


Materials 

Planes 

2/7,  deg 

IJ.  cm  ' 

(5,  10  ’  in,  for  absorption  of: 

95% 

993  0 

Al 

(311) 

139.2 

402.3 

1.37 

2.11 

Ti 

a  -  hep 

(201) 

136.3 

2737.6 

0.200 

0.296 
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Figure  13.  d-spacing  for  (200)  plane,  7475  Al,  using  CuKa,  at  (lr  =  0,  as  a  function  of 
distance  from  the  centerline. 
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Figure  14.  d-spacing  for  (400)  plane,  7475  Al,  using  CuKa,  at  i^  =  0,  as  a  function  of 
distance  from  the  centerline. 
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4.5  d-spacing  calculations 


The  strain  distribution  illustrated  in  Figure  8  was  used  to  obtain  the  d-spacing  variation  as  a  func¬ 
tion  of  the  distance  from  the  centerline.  This  was  done  for  a  7475  AI  specimen  using  CuK«  radi¬ 
ation.  This  profile  was  used  primarily  to  simulate  the  line  profile  on  the  flat  surface  (A)  of  the 
notched  specimen  Details  regarding  this  simulation  are  described  in  this  chapter.  As  the  x-ray  re¬ 
sponse  in  terms  of  the  line  profiles  was  studied  on  the  flat  surface  of  the  specimen,  the  d-spacing 
profiles  were  also  obtained  for  the  same  surface. 

Consider  an  x-ray  beam  incident  on  surface  A  of  the  specimen  as  shown  in  Figure  3.  The  inter- 
planar  spacing  d  of  the  diffracting  planes  parallel  to  the  surface,  result  from  the  strain  normal  to  the 
'ur^ace.  Ihis  strain  for  ip  =  0  is  the  strain  in  the  i.  direction,  i  c.,  c,  or  r.,,.  Ihc  d-spacing  profiles 
for  three  orders  of  reflections  namely  (200),  (400)  and  (333),  resulting  from  c,  strain  arc  illustrated 
in  Figures  13-15.  The  value  of  d,  for  these  reflections,  were  obtained  from  Pearson's  Handbook  (26|. 
These  profiles  were  obtained  from  the  strain  distribution  of  Figure  8  and  equation  (4.02).  Figure 
13  exhibits  a  gradient  which  is  steep  in  the  region  close  to  the  notch  and  decreases  gradually  to  zero 
at  the  largest  distance  away  from  the  notch.  This  can  be  attributed  to  the  direct  proportionality  of 
the  d-spacing  with  the  strain,  which  displays  a  similar  behavior  (I'igure  8) 
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4.6  Simulation  of  intensity  profiles 


4.6.1  Introduction 


The  most  common  kind  of  data  one  obtains  from  a  diffractometer  is  a  graph  showing  the  variation 
of  the  diffracted  intensity  as  a  function  of  W.  Ilieorctically  the  diffracted  intensity  is  defined  to  be 
a  step  function.  However  in  reality,  factors  like  the  presence  of  a  mosaic  structure  (small  regions 
slightly  misoriented  from  one  another)  in  the  crsstal,  size  of  the  crystal,  presence  of  nonuniform 
elastic  stresses  and  strains  and  instrumental  effects  may  lead  to  a  broadening  of  the  diffracted  in¬ 
tensify.  Such  a  broadened  diffracted  intensity  is  known  as  a  line  profile. 

This  section  describes  the  simulation  of  intensity  or  line  profiles  resulting  from  the  d-spacing  pro¬ 
files  illustrated  in  Figures  13,  14  and  IS.  The  simulation  was  carried  out  for  an  x-ray  beam  incident 
on  surface  A  of  the  notched  specimen  (Figure  3).  Intensity  profiles  were  simulated  for  diffraction 
from  each  grain  in  the  region  between  the  notches,  for  a  beam  of  finite  size.  The  profiles  were 
simulated  to  study  stress  gradient  effects,  and  thus,  the  effects  of  a  non  uniform  macrostrain  on 
stjrface  A  of  the  specimen.  The  same  simulation  was  not  done  for  surface  B  at  the  bottom  of  the 
notch.  This  is  because,  for  surface  B.  for  a  finite  beam  size  used  in  the  simulation  (and  also  in  the 
experimental  work),  the  strain  gradient  relative  to  the  depth  of  penetration  of  the  x-ray  beam  (dis¬ 
cussed  in  section  4.3)  is  negligible.  Hence,  the  effect  of  this  gradient  on  the  line  profile  would  be 
insignificant.  All  the  results  presented  in  this  section  were  obtained  by  running  the  simulation  for 
a  7475  A1  notched  tensilr  specimen. 
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4.6.2  Theory 


The  d-spacing  profiles  from  figures  13-15  were  used  in  the  simulation  of  the  intensity  profiles.  The 
intensity  profiles  corresponding  to  the  d-spacing  variations  were  generated  using  a  computer  routine 
written  specifically  for  this  purpose.  It  should  be  noted  that  the  strain  distribution  used  as  the  basis 
for  these  calculations  corresponds  to  y  =  0.0,  ic.,  along  the  centerline  of  the  notches.  The  variation 
in  the  stress  distribution  between  y  =  0.0  and  0,1  has  already  been  shown  to  be  small  (Figure  7). 

The  20  values  corresponding  to  the  d-spacings  calculated  in  the  last  section  can  be  comr"^!cd  from 
equation  (2  01)  as 


0  =  sin  ' 


(4.22) 


The  continuous  variation  of  the  d-spacing  pves  rise  to  an  intensity  profile  spread  over  a  20  range 
according  to  Bragg's  law.  However,  the  simulated  profile  includes  broadening  due  to  the  sample 
and  the  instrument.  This  can  be  accounted  for  by  the  use  of  a  Pearson-Vll  |30]  function  This 
function  can  be  used  widely  for  fitting  profiles  of  Cauchy  and  Gaussian  shapes,  including  the 
modified  Ixirentzian  shape.  The  P-VII,  as  it  is  usually  called,  is  given  by, 


>  = 


1  -t- 


(20  -  20^)^ 


ma 


(4.2.1) 


where,  the  parameters  Y,  and  20,  indicate  the  intensity  and  the  position  of  the  peak.  The  integral 
of  (3')  with  infinite  limits  with  respect  to  0  is  noimaiir.cd  if 


r(m) 

<’  r(m-I/2) 


(4.24) 


where,  a'  and  m'  are  two  adjustable  parameters  that  describe  the  peak  shape  and  the  width  re¬ 
spectively.  The  full  width  at  half  maximum  (fTVHM)  can  be  shown  to  he 
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W=2a[m{2'l’^  -\)]'<^  (4.25) 

When  m  =  1, 2  and  oo,  the  P-VIl  becomes  a  Cauchy,  modified  lorcntzian  and  Gaussian  cunc 
respectively. 

In  order  to  obtain  a  smooth  variation  in  the  intensity  profile,  the  d-spacing  profile  calculated  from 
the  strain  distribution  was  expressed  as  a  series  of  five  polynomials  This  five  segment  approxi¬ 
mation  of  the  d-spacing  profile  is  shown  in  Figure  16.  The  five  segment  approximation  corre¬ 
sponding  to  the  higher  order  reflections  arc  also  shown  in  Figures  1 7  and  1 8. 

4.6.3  Intensity  profile  calculations 


A  computer  simulation  approach  to  obtain  composition  profiles  has  been  done  earlier  in  diffusion 
studies  by  Ilouska  and  co-workers  This  has  been  reviewed  in  an  article  by  llouska  1.51).  In  the 
approach  by  Houska,  the  variation  of  the  composition  and  depth  was  considered  However,  in  the 
present  model,  as  the  d-spacing  variation  with  depth  is  negligible,  the  d-spacing  variation  along  the 
x-axis,  as  opposed  to  depth,  was  considered 

A  computer  program  was  written  to  calculate  the  intensity  profiles.  The  program  used  as  input  the 
d-spacing  of  the  planes  in  the  unstressed  state  (rf.)  for  the  specific  reflections  considered  which  in 
the  present  case  were  (200),  (400)  and  (3,13).  ITic  height  above  the  centerline  of  the  notches,  y,  at 
which  the  stress  distribution  had  to  be  determined  was  also  input  into  the  program.  The  program 
is  capable  of  determining  the  stress  distribution  and  the  line  profiles  for  different  values  of  y. 
However,  the  intensity  profiles  have  been  illustrated  for  a  y  value  of  0.0  only.  The  other  inputs 
correspond  to  the  coefTicicnts  of  the  polynomial  equations  used  in  the  five  segment  approximation 
of  the  d-spacing  profile.  Sample  input  files  are  shown  iii  Appendix  R. 
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Figure  16.  Five  segment  approx  of  d-spadng  for  (200),  7475  Al,  CuKa,  at  =  0  as  a 
function  of  distance  from  the  centerline 
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Figure  17, 


Five  segment  approx  of  d-spacing  for  (400),  7475  AJ,  CuKa,  at  \l> 
function  of  distance  from  the  centerline 


0  as  a 


Distonce  Iron  Ih*  centerline  (inches) 


Figure  18.  Five  segment  approx  of  d-spacing  for  (333),  7475  Al,  CuKa,  at  \l/ 
function  of  distance  from  the  centerline 
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The  program  is  divided  into  the  following  four  sections  ; 

•  The  main  program  calculates  the  intensity  profile. 

•  Subroutine  STRAIN  calculates  the  stress  and  strain  values  using  the  equations  presented  in 
chapter  3. 

•  Subroutine  DX2TT  calculates  the  d-spacing  variation 

•  Subroutine  DVSSIN  calculates  the  d  vs  sin*  (//values  for  the  model  developed  earlier. 

The  main  program  first  divides  the  representative  area,  which  in  this  case  is  the  region  along  the 
centerline  of  the  notches,  into  many  fine  strips  along  the  x-axis.  As  the  d-spacing  gradient  in  each 
strip  is  of  a  small  magnitude,  the  d-spacing  in  each  fine  strip  is  approximated  as  a  constant.  The 
intensity  for  each  strip  was  first  calculated.  Ibis  intensity  was  broadened  using  the  Pearson  N  il 
function  (equation  4.23).  The  broadened  intensity  was  then  superimposed,  by  means  of  a  sum¬ 
mation  loop,  to  obtain  the  simulated  intensity  profile.  This  summation  is  a  summation  of  intensi¬ 
ties.  For  the  sort  of  summation  done  in  the  present  situation,  the  coherence  length  of  the  subgrain 
was  assumed  to  be  small.  The  depth  of  penetration,  d-spacing  gradient  and  the  subgrain  size  of  the 
material  are  factors  which  affect  the  calculations.  The  depth  of  penetration  of  a  CuK<x  radiation  in 
A1  has  already  been  found  to  be  small  (Table  4).  For  this  depth  of  penetration,  the  variation  in  the 
d-spacing  due  to  the  .strain  c„  is  small.  These  two  factors  and  the  as.sumption  of  a  small  subgrain 
size  help  in  justifying  the  summation  of  intensities  in  our  calculations.  The  intensity  profile  based 
on  the  above  calculation  is  .shown  in  Figure  I'f  I  hc  program  was  then  run  for  higher  order  re¬ 
flections.  namely,  (400)  and  (.3.3.3).  Ibe  profiles  generated  for  these  arc  also  illustrated  (Figures 
20-21). 

Broadening  observed  in  any  intensity  profile  is  actually  a  convolution  of  four  different  broadening 
functions  Ibese  are  :  broadening  due  to  non-uniform  strain,  dislocation  strain  broadening,  particle 
size  broadening  and  instrumental  broadening  Ibc  intensity  profiles  illustrated  in  Figures  (19-21) 
depict  a  clear  asymmetry.  In  order  to  know  the  contribution  of  each  of  the  broadening  function 
towards  the  asvmmclry.  a  line  shape  analysis  of  the  profiles  is  required  As  the  method  by  which 
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the  intensity  profile  was  generated,  is  a  linear  summation  over  a  scries  of  fine  strips  along  the  x-axis, 
information  of  the  kind  required  in  the  present  situation  can  be  obtained,  by  plotting  the  results 
due  to  the  summation  over  a  smaller  range.  This  result  is  shown  in  Figure  22.  In  this  case,  half  the 
region  between  the  notches  was  divided  into  approximately  1200  linear  elements  and  the  intensity 
due  to  each  element  was  obtained  from  the  computer  program.  liach  peak  was  then  broadened 
using  the  Pearson  VII  function  discussed  earlier  in  this  chapter.  The  linear  elements  were  then 
summed  twice.  The  results  were  first  summed  over  the  entire  region  between  the  notches  (Figure 
3),  to  generate  the  profile  revealing  the  asymmetry  in  Figures  19,  20  and  21 .  Fhc  second  summation 
was  over  separate  regions,  which  covered  one  tenth  the  region  between  the  notches  along  the  notch 
centerline  on  the  flat  surface  (A),  each  of  which  is  equivalent  to  the  width  of  the  slit,  n.5mm  (0.02  ). 
used  in  the  TFC'  x-ray  system  leioking  at  the  profiles  generated  by  these  ten  segments  (Figure  22). 
it  is  clear  that  the  profile  with  the  maximum  peak  corresponds  to  the  region  farthest  away  from  the 
notch,  and  thus,  a  low  strain  gradient  As  each  peak  corresponds  to  a  normalised  distribution 
function,  the  area  under  the  curve  is  taken  to  be  a  constant. 


4.6.4  Discussion  of  the  intensity  profiles 


The  intensity  profiles  were  found  to  be  very  sensitive  to  variations  in  the  d-spacings,  ic.,  to  the  d- 
spacing  gradient.  Small  variations  in  the  inter-planar  spacings  lead  to  large  changes  in  the  slope, 
which  give  rise  to  large  changes  in  the  shape  of  the  intensity  profile  With  reference  to  Figure  19, 
the  peak  of  the  intensity  profile  is  found  to  correspond  to  a  small  d-spacing  gradient.  This  is  be¬ 
cause,  for  low  gradients,  the  x-ray  beam  covers  a  larger  region  of  material  and  thus  a  greater  number 
of  grains.  It  should  be  noted  that  the  profiles  simulated  here  correspond  to  a  d-spacing  variation 
along  the  x-axis,  on  surface  A  (f'igurc  3)  of  the  specimen.  Hie  d-spacing  at  each  position  along  the 
x-axis  is  a  result  of  the  strain  in  the  7  direction  and  thus  for  =  0  Hence,  for  this  simulation,  the 
absorption  for  each  element  along  the  x-axis  of  the  x-rays  incident  on  surface  A  was  assumed  to 
be  constant. 
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Figure  19.  Intensity  profile  for  (200)  plane.  7475  Al,  using  CuKa  radiation  through  a 
rectangular  slit. 
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Figure  20- 


Intensity  profile  for  (400)  plane.  7475  Al.  usmg  CuKa  radiation  through  a 
rectangular  slit. 
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Figure  21.  Intensity  profile  for  (3331  plane,  7475  A1 ,  using  CuKa  radiation 
through  a  rectangular  slit. 
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Figure  22.  Intensity  profile  for  (200)  plane,  7475  Al,  using  CuKa,  from  ten  segments,  each 
segment  of  width  0.02'. 
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In  the  calculation  of  the  intensity  profiles,  as  stated  earlier,  the  region  in  between  the  notches  on 
surface  A  was  divided  into  a  number  of  strips.  The  si/e  of  the  strip  used  in  the  computation  was 
small  compared  to  the  grain  size  of  the  material  The  depth  of  penetration  of  the  x-ray  beam  was 
found  to  be  small  too  (refer  Table  4).  The  coherence  length  of  a  subgrain  which  in  this  case  is  the 
column  length,  is  assumed  to  be  small  for  this  depth  of  penetration.  The  profiles  shown  in  Figures 
19-21  were  obtained  by  summing  up  the  intensities  from  each  column.  The  asymmetry  in  the  pro¬ 
file  is  best  understood  by  treating  each  domain  in  the  region  between  the  notches  on  surface  A  in 
terms  of  a  column  of  cells 

Consider  a  column  of  cells  as  illustrated  in  Figure  2.^  wherein,  ai,  a;  and  ai  arc  the  directions  along 
the  cel!  edges  F.ach  domain  can  be  expressed  in  terms  of  a  column  of  cells  along  the  aj  direction, 
perpendicular  to  the  reflecting  planes,  i  and  i'  represent  different  columns.  From  the  expression  for 
integrated  intensity  |32],  it  is  seen  that  contributions  to  the  intensity  arc  only  from  pairs  of  cells  for 
which  =  m,  and  zwj'  =  m2.  This  implies  that  there  is  no  contribution  to  the  integrated  in¬ 
tensify  in  each  column  due  to  cross-over  terms.  In  other  words,  cells  from  different  columns  do  not 
contribute  to  the  integrated  intensity  in  each  column.  As  the  x-ray  beam  not  only  scans  all  columns 
in  a  single  grain  or  subgrain,  all  other  grains  having  similar  orientations  also  contribute  towards  the 
intensity  profile.  The  intensity  profiles  obtained  by  the  simulation  ha\c  been  by  a  superposition 
of  the  intensities  from  different  columns  and  subsequently  from  different  grains. 

1  his  theory  for  the  calculation  of  the  intensity  relating  to  a  column  of  cells  is  for  a  characteristic 
d-spacing.  Ihc  program  calculates  the  intensity  for  a  group  of  cells  in  the  domain  corresponding 
to  a  particular  d-spacing  Different  groups  of  cells  over  the  domain  Icail  to  a  variation  in  the  d- 
spacing  The  overall  profile  is  obtained  by  a  summation  over  all  the  groups  of  cells  in  the  entire 
domain. 

The  profiles  corresponding  to  the  lower  order  reflections  depict  a  very  narrow  2/?  range.  The  profile 
for  the  (.33.3)  reflection  shows  a  2d  range  of  approximately  2'.  The  variation  in  this  profile  is  con¬ 
siderably  sharp,  with  the  predominant  contribution  to  the  asymmetry  being  from  the  tail  portion. 
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The  effects  of  the  intensity  profiles  may  not  he  easily  measurable  on  the  1  PC'  x-ray  system  being 
used  in  the  experimental  work.  They  can  however  be  measured  on  an  x-ray  diffractometer  with  a 
higher  resolution  using  a  Ka,  monochromator. 

Consider  the  simulated  profiles  illustrated  in  F'igurcs  (19-21).  These  profiles  are  plotted  for  three 
different  orders  of  reflection  (hkl).  For  higher  orders  of  reflection,  the  inter-planar  spacing  decreases 
according  to  the  equation, 

d=-  ^-=1  (4  26) 

where  a  is  the  lattice  parameter  of  the  material  from  which  diffraction  occurs.  Higher  order  of  re¬ 
flections  correspond  to  a  lower  d  value  and  hence  a  larger  20.  as,  the  inter-planar  spacing  and  P  have 
an  inverse  dependence  (Bragg's  Law). 

The  peak  shape  depends  on  the  size  of  the  source,  finite  divergences  of  the  slits  at  the  x-ray  tube 
and  the  receiving  slits,  and  the  fact  that,  in  a  diffractometer,  a  flat  specimen  is  tangent  to  the 
focussing  circle  at  only  one  point.  All  x-ray  measurements  on  the  I  fT!  x-ray  system  are  on  the  higli 
angle  side.  At  high  angles,  all  the  above  factors  become  less  important.  Also  at  high  angles,  the  re¬ 
solution  in  Afl  is  larger  compared  to  the  resolution  at  low  angles  This  is  evident  from  equation 
2.02b  which  is 


Afl  =  — tan  0  =  —  r.  tan  P 
dhki 

Tbese  high  angle  measurements  are  confined  to  a  2/1  range  of  140"  -  170".  The  profile  for  (.13.7) 
reflection  is  for  a  high  angle  diffraction,  i.e.,  20  ~  160°.  This  profile  shows  an  asymmetry  over  a 
20  range  of  2'.  The  profiles  for  (200)  and  (400)  are  for  diffraction  at  a  2/1  of  ~  45"  and  90"  re¬ 
spectively. 
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In  an  actual  experimental  measurement,  a  broadened  intensity  profile  of  a  form  similar  to  I'igures 
(19-2!)  is  obtained.  The  broadening  due  to  the  instrument  can  be  obtained  from  measurements  on 
a  well  annealed  sample.  Measurements  made  on  a  sample  with  very  few  defects,  using  a  system 
whose  instrumental  broadening  is  known,  yields  the  sample  or  particle  size  broadening.  A  measured 
profile  is  usually  from  a  deformed  sample  in  a  real  instrument.  This  profile  is  a  convolution  of  the 
the  broadenings  due  to  the  sample  and  the  instrument,  and  also  the  effects  due  to  a  non-uniform 
macrostrain.  The  broadening  due  to  the  non-uniform  macrostrain  can  be  determined,  if  the  factors 
due  to  the  sample  and  the  instrument  are  known. 

From  the  experiments  done  on  the  TIT'  x-ray  system,  the  full  Width  at  Half  Max  (FAVMMi  for 
diffraction  on  Aluminum  specimens  was  ~  2''  -  T.  'llic  FWIIM  is  really  a  measure  of  the  resolution 
of  the  instrument.  (Considering  the  fact  that  one  of  the  effects  exhibited  by  the  simulation  is  of  the 
same  order  ir  2'  the  effect  may  be  measurable  experimentally,  but  with  some  difficulty.  This  effect 
should  be  easily  measurable  on  a  diffractometer  with  a  higher  resolution 
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5.0  Conclusions  and  Recommendations  for  Future  Work 


Based  on  the  results  of  the  analvtical  model,  the  following  conclusions  can  b-e  drawn 

•  Calculations  related  to  the  stress  distribution  model  for  an  infinite  specimen  can  be  extended 
to  a  specimen  of  finite  gcomctr>'.  The  error  encountered  due  to  this  approximation  is  of  the 
order  of  3  percent. 

•  The  calculations  regarding  the  variation  of  the  inter-planar  spacing  with  the  i//  angle,  illustratcil 
as  a  plot  of  d  vs  sinV  .  can  be  used  as  a  basis  for  experimental  measurements.  The  variation 
in  the  intcrplanar  spacing  arc  easily  measurable  experimentally  The  tabulated  values,  calcu¬ 
lated  for  a  point  at  the  notch  tip,  can  be  verified  experimentally.  In  this  context,  attaining 
precision  in  the  x-y  positioning  of  the  x-ray  beam  on  the  specimen,  is  of  prime  importance. 

•  A  study  of  the  effect  ,of  the  depth  of  penetration  of  the  x-ray  beam  on  the  stress  gradient 
occuring  on  the  surface  reveals  that  residual  stress  measurements  arc  feasible  along  the  edge 
of  the  specimen 
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•  The  intensity  profiles  simulated  arc  due  to  the  strain  gradient  on  the  surface  They  arc  a  result 
of  a  non-unifonn  macrostrain  in  a  notched  specimen  subjected  to  tensile  loading. 
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Appendix 

Effect  of  the  radius  of  curvature  of  the  sample  ou  the  error 
in  d-spacings 


For  a  sample  with  a  radius  of  cutA  alurc  the  total  shift  in  the  centroid  of  the  diffraction  peak 
has  been  found  to  he 


2  a^s\i?e 


1  + 


rgc(l 

2r5  sin(0  4-  i^) 


(A  -  1) 


where  rg^,  is  the  radius  of  the  goniometer  circle  with  the  center  on  the  point  where  the  x-ra\  beam 
strikes  the  sample  and  passes  through,  the  x-ray  tube  focal  spot,  rpj.p^  is  the  distance  between  the 
centers  of  the  goniometer  and  the  detector,  the  detector  being  parallel  to  a  tangent  of  the  goniometer 
circle,  6  is  half  the  Bragg  angle,  r^  is  the  radius  of  curvature  of  the  sample,  and  n  is  half  the  angular 
beam  divergence.  Note  that  convex  samples  have  r^  >  0  and  concave  samples  have  r^  <  0.  The 
angular  beam  divergence  for  a  point  source  and  a  source  of  finite  height  is  shown  in  Figures  A-I(al 
and  A- 1(b)  respecliscly,  in  which  d  is  the  width  of  the  slit,  u  is  the  distance  between  the  two 
apertures,  v  is  the  distance  of  the  exit  pinhole  from  the  source,  h  is  the  height  of  the  source,  and 
2'»|  and  2^2  maximum  angles  of  divergence  for  the  two  cases 
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I  he  source  pertaining  to  the  present  investigation  is  square  in  shape  with  h  appro ximatelv 
being  equal  to  2  mm  (Figure  A- 1(b)).  The  divergence  in  radians  is  then  given  by, 

_  (d/2)  _  d 
“2  (u/2)  u 


Hence, 


(A -2) 


Tsing  equation  (.A-1)  and  difTcrentiating  equation  (2.01)  and  solving,  the  error  in  the  d-spacing  is. 


d 


r  ')  1 

2  a^cos^O 


i  + 


‘gc' 


-r) 


2r,.  sin(fl  +  f 


(A  -  .^) 


Fquation  shows  that  the  shift  in  d-spacing  involves  a  term  which  is  inversely  proportional  to 
the  radius  of  curv  ature  of  the  sample.  Substituting  the  following  knowm  values  for  the  diffractometer 
being  used  in  the  experimental  work,  i.e., 
d  =  0.5mm, 
u  =  5()mm. 

a  =  0.(1  i  radians  =  0,6^, 
rg^  =  2n0mm  and 

'pspc 

Fq.  (A-.f)  becomes. 


d 


=  6  67  X 


in 


-5 


cos  d 


sm(0  -t- 1/') 


1 


98.01 

rj  sin(0  +  4') 


(A -4) 


The  values  of  in  parts  per  million  for  different  i  angles  corresponding  to  each  of  the  20 

bracket  used  arc  shown  in  Table  A- 1  f-quations  (.A-1)  and  (.A-4)  presented  here  arc  denvations  by 
Hendricks  [  t,t| 
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Tabic  A-1.  Estimate  of  the  error  in  d-spacing  values  for  differing  radius  of  curvature. 


Radius  of 
curvature 

Material 

V'  =  0° 

10° 

Ad  ^  \ 

-j-  (ppm) 

20°  30° 

45° 

3. 175mm 

Al' 

-62.9 

-60.0 

-62.9 

-72.5 

-109.9 

li'' 

-250.1 

-219.1 

-211.2 

-223.3 

-291.7 

6,35mm 

Al 

-30.4 

-29.0 

-.304 

-35.1 

-5.3.5 

li 

-121  1 

-106  0 

-102.0 

-lOR  0 

-141.5 

■  20  =  160’  and  =  =  142’ 

I  hcsc  values  have  been  calculated  for  two  different  radii  of  curvature,  r^  =  -t  mm  and  r^  =  -b 
mm  respectively,  the  negative  signs  implying  a  concave  sample  surface.  The  error  in  the  residual 
stresses  measured  corresponding  to  the  above  error  in  d-spacing  can  also  he  estimated  From 
equation  (2.06''  "c  have, 


Ad  (1  +  ’) 


■  2 


=  — p - sin  i/'  -  —  (o  I  +  02)  (A  -  5) 


[differentiating  both  sides  of  the  above  equation  with  respect  to  sin^i/'  ,  vve  have. 


d(  sin^i//) 


I  Ad  ,  Ad 

^  d  ’'I' 2  ^  d  ’'I']  ^  I  +  V 

•  2  -2  f- 

sin  «/'  2  —  sin  i/'  I 


(A -6) 


In  I  quation  (A-b),  the  values  of  the  change  in  stress  da^  can  be  calculated  knowing  the  values  of 

from  Table  A-1.  These  calculations,  done  for  both  7475- A1  T651  and  ri-bAl-4V  allovs.  arc 
d 

show  n  in  I  able  A-2  The  results  of  T  able  A-2  indicate  that  there  is  no  significant  change  in  the 
stress  values  for  both  the  sharp  and  blunt  notches,  l-.quations  ^.^-l)  and  (A-4)  arc  consenative 
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overestimates  of  the  effect  dealt  with  in  this  section.  Therefore,  the  results  of  T  ables  A-1  and  A-2 
are  overestimates  too  [.T4]. 

Table  A-2.  Estimate  of  the  error  in  stress  values  for  differing  radius  of  curvature. 

Material  Aai 

(in  l^i) 

A1  0,<?5 

Ti  2.50 

A1  0.47 

Ti  1.2.7 
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Listing  of  program  INTENSITY 


C  INTENSITY 

C 

(;***,**************************************************,**»**«««*.*.,«« 

C  THIS  IS  A  PROGRAM  TOR  INTENSITY  BAND  SIMULATION  FROM  A  NOTCH  SPF  CIMEN 
C  or  7475AL-T651,  FOR  AN  X-RAY  BEAM  IRRADIATING  I  Hi:  SURI  ACE  OE  THE 
C  SPECIMEN  IN  THE  REGION  BETWEEN  THE  NOTCHES. 

C 

C  INPUT  PARAMETERS 
C 

C  DO  IS  HIE,  D-SPACING  01  THE  UNSI  RESSED  PI  AM  S 
C  ALAMDA  IS  THE  WAVELENGTH 
C  YI  IS  HIE.  HEIGHT  ABOVE  THE  BASE:  OI  THE  NOI  CII 
C  X2  IS  THE  POSITION  ALONG  THE  X-AXIS 

C  A(I),  B(I)  AND  C(I)  ARE  THE  VALUI-S  OE  HIE:  COE.EEICIIMS  IN  HIE 
C  POLYNOMIAL  USED  TO  EXPRESS  HIE  D-SPACING  PRO!  II  E 
C  SX.  SY  ARE  THE  STRESSES  IN  THE  X  AND  Y  DIRECTION  RI  SIM  C  l  IVl  I  Y 
C  EX,  EY,  AND  EZ  ARE  THE  CORRESPONDING  STRAINS 
C  IT  IS  the:  theta  value  CORRESPONDING  EO  THE  D-SI’ACING 
C  NX  IS  THE  NUMBER  OF  POINTS  AI  WHICH  SI  RAIN  RESUI  IS  ARE  OB  EAINI  D 
C  NT  IS  THE  1  OTAL  NUMBER  OE  POIN  IS  GENERA  l  ED 
C  SPSI  IS  THE  SIN**2  VALUE 

C  DOS,  DOE.  D90S,  AND  D90E  ARE  THE  D-SPACING  RESUL  TS  ON  SUREACE  (S) 

C  AND  EDGE  (E)  RESPECTIVELY  FOR  PHI  =  0  DEG.  AN  DEG.  RESPECHVELY 
C  START  IS  THE  TWO-THETA  VALUE  (IN  DEGREES)  OE  THE,  LOW  ANGLE  SIDE 
C  or  HIE  INTENSITY  BAND  TO  BE  SIMULATED 
C  DELTT  IS  THE  INCREMENT  IN  IWO  HIE:TA. 

(^********«********«*****t«***************«»*********«**«****« »**«**» 

C 

IMPLICIT  REALM(A-H,0-Z) 

DIMENSION  D(2000),X(2000),7T(2000),X2(  100) 

DIMENSION  P(2000),UBP(2000).SINrEN(2000),AINTEN(2000),TIIErA(200n) 

DIMENSION  EX(100),EY(100),EZ(IOO),SX(100),SY(100) 

DIMENSION  DOS(50),D90S(50).DOE(50),D90E(50),SPSI(50) 

CALL  ERRSET(208.256,-1) 

ALAMDA  =  1,54178 
PI  =  ACOS(-l.O) 

RAD=  2.0‘180.0/PI 
C 

C  THi:  F  OLLOWING  ARE  HIE  READ  SrA  EE:MI'NIS  IN  HIE,  PROGRAM 

c: 

READ(5.‘)D0 

rf:ad(5,‘)yi 

rf:ad(5,*)ai,bi,ci 

READ(5,*)A2,B2,C2 

READ(5,*)A.1.B3,C.'' 

RI:AD(5,‘)A4,B4,C4 

RE:An(5,‘)A5,B5,C5 

C 
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C  THF  FOl.LOWING  ARF  THF  FORMA  I  S  I  ATFMFN  FS  IN  I  III'  PROGRAM 

C 

1  FORMAT('X',8X,'SX’.8X,'SY',') 

2  FORMAT(/5X,'X2',8X.'EX',8X.'EY'.8X/FZ',/) 

3  FORMATS  D',8X,'X',8X.TWO-Tl!ETA  ,/) 

4  FORMAT(/TWO-THETA'.8X.'INTENSn  Y  ,/) 

5  FORMAT(2X.F10.4,2fF10.3)) 

6  FORMAT(2X,F10.4,3(F10.3)) 

7  rORMAT(2X,F12.6.F124,F13.6) 

8  FORMAT(2X.F13.6,E13.6) 

9  rORMAT(2X,F5.2.4(F13.5).5X,F5.2) 

11  FORMAT(3X,'SPSI  ,5X.'D0SPSl’,8X.'D0FPSr,7X,'n9nSPSr.5X,'D901PS 
*r,7X/Y',/) 

Q^t********************************************************************* 

c 

CALI  S1RAIN(Y1,X2,FX.FY,EZ.SX.SY,NX) 

CAI  .1 .  DX2n  (D0,NX.I-Z.D,T  I ) 

CAFF  nVSSlN{D0,X2,FX,FY.I'Z.SPSFn0S.DnF„D9nS,l)9ni  1 
NP=  1270 
X(i)  =  n.n 

c 

^-'«««*««««***«**««***t*«*««*»4***««***«**««*«*t************i«*«*«***i«««< 

C  THF  FOFl  OWING  IX)  FOOP  GIVFS  Fill  SI  T  OI'  FgFA  I  IONS  FOR  I  III  O-SPACING 
C  CURVE  FI  I 


DO  10  1=  FNT 

IF(X(1).FE.0,04)  D(I)  =  A1*(X(I)*»2)+  BI*X(I)+  Cl 
IF(X(I).GT.0.04  .AND,  X(I).FE.0.08)  D(I)  =  A2*(X(1)**2)  +  B2‘X(1)  +  C2 
IF(X(1).GT.0,08  .AND.  X{I).FE.0.1S)  D(l)  =  A3*(X(I)‘»2)  +  B3*X(I)  +  C3 
IF(X(I).GT.0.15  .AND.  X(I).EE.0.21)  D(I)  =  A4*(X(I)**2)  +  B4*X(1)  +  C4 
IF(X(I).GT.0.21  AND.  X(I).LE.0.26)  D(I)  =  A5*(X(I)»*2)  +  B.MX(I)  +  C5 
1T(I)  =  ASIN(AFAMDA/(D(1)‘2.0))»RAD 
X(I+  l)  =  X(l)  +  2.nF-04 
in  CONTINUE 
C' 

C  FlIF  FOFFOWING  DO  L  OOP  CORRESPONDING  FO  STA  FI.MF  N  I  #  100  CJI  NL  RA  FI  S  I  HI 
C  IN  FFNSn  Y  DUE  FO  EACH  OF  THE  FFN  SI  GMl  NFS  IN  Fill  RI  (HON  BI  I  WI  I  N  Fill 
C  NOTCIIF  S 

r 

N  =  n 

DO  inn  1=  1.10 
M  =  I  +  N 
N  =  N+  127 
C 

C  FOOP  20  CAI  CUI  A7FS  FUF  IN  I  FNSIFV  FOR  FAC'II  FAYI  R  IN  A  FI  A  I  SAMPI  I 
C 

IX)  20  1=  M,N 

20  PH)  =  (X(1+ l)-X(I))/0.25 

(.• 

C  INITIAFISING  Till-  F’NBROADFNFD  IN  FFNSFI  V  AND  THE  IN  FI  NSII  V  SIMUI  A  l  l  D 
C  AT  THF  1  A(l)  1  HF  I  A(l)  IS  Fill  1  WO- 1  HI  FA  POSI  FION  A  I  WHK  H  FHI 
C  INFI  NSIFII.S  ARF  TO  Bl  SIMUI  All  D  SINFFN  IS  FHF  INFI  NSFFV  DUt  FO 
C  EACH  SFGMFNF 
( 
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DO30I  =  M,N  ! 

i]BP(i)  =  o.n 
30  SINTF=N(l)  =  0.0 
C 

AN  =  FLOAT(N) 

DELTT  =  0.3/127.0 

START  =  TT(N)  -  0.10  ! 

C 

DO  40  I  =  M,N 

40  THETA(I)  =  START +DELTr*(l-M) 

THE  FOLLOWING  LOOP  CALCULAl  l-S  1 1  IF.  IN MLNSH  Y  FOR  EACH  POINT  OF  THE. 
INTENSITY  BAND  BEING  SIMULATED. 

DO  60  K  =  M,N 

IF(K.EO  1  OR.K.E0.NP)GO  TO  108 
AM  =  ( rHETA(K)  +  THETA(K-l))/2  0 
ULI  =  (THFTA(K)  +  THETA(K+  1))'2.0 

GW  AND  GM  ARE  TWO  ADJUSTABI  i:  PARAMF.I  I  RS  DI  NO  I ING  l  llF  WIDl  II  OF  HIE 
GUASSIAN  PEAK  AND  THE  POWER  FOR  HU .  GLASSIAN  DISlRIBrilON  RESPl  CUVEL' 

108  GW  =  0.01 
GM  =  2.5 


THIS  LOOP  BROADENS  HIE  CALCULA  Ii;!)  INI  ENSM  Y  OI  I.ACII  LAYl.K  USING  HII- 
PEARSON  VII  FUNCTION.  THESE  ARE  I  I11:N  SUPI  RIMPOSI  D  K)  GI  1  HIE  SIMULATED 
INTENSITY. 

DO  50  J  =  M,N 
TTJ  =  (TT(J)  +  TT(J  +  1  ))/2.0 
DIFF  =  THETA(K)-TTJ 

YNORM  =  (l.0/(SQRT(GM*PI)*GW))*(GAMMA(GM)/GAMMA(GM-0..8)) 

BR  =  YNORM/(1.0  +  (DIFF**2)/(GM*(GW»*2)))»*GM 
IF(K.EQ.  l  .OR.K  EQ.  NP)GO  TO  50 

IF(  rr(J)  GE.ALI  -AND.  TT(J)  .LT.ULI)UBP(K)  =  UBP(K)  +  P(.I) 

)0  SINTEN(K)  =  SINTEN(K)  +  P(J)*BR 
>0  CONHNUF 

THIS  WRITE  LOOP  WRITES  THE  IN  I  ENSU  Y  VALUES  DUE  3  0  EACH  OF  THE  TEN 
SEGMEN  I  S  ALONG  WITH  THEIR  CORRESPONDING  TWO-HIF  I  A  POSITION 

WRITE(6,4) 

DO  45  1=  M,N 
SINrF.N(I)  =  SINIFN(I) 

WRH  F,(6,8)THr  rA(I).SINTEN(I) 

15  CONTINUE 
[)0  CONTINUE 

»**t*****«*********««********************«*****************t««****** 


THIS  SECHON  USES  THE  SAME  PRINCIPLES  AS  THE  PREVIOUS  SI  ('HON 
IN  THE  CAI  CULATION  OF  IN  I  ENSU  Y.  HIE  ONLY  DILI  ERI  NCF  IS  THAT 
HERE  Tin  OVERAL!  INI  ENSU  Y  OVf  R  HIE  ENHRE.  DOMAIN  IS  CAI  (  LI  AH  I) 
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C  THIS  OVr.RAI  1  INTFNSITY  IS  THRMPn  AS  AIM TN 
C' 

N  =  NR 
M=  I 

DO  120  1  =  M,N 
120  P(I)  =  (X(1+ l)-X(I))/0.25 

C 

DO  130  1  =  M,N 
UBP(l)  =  0.0 
130  A1NTHN(I)  =  0.0 
C 

ANP  =  ri,OAT(NP) 

DE1TT  =  0.3/ANP 
STAR  !  =  1  l(N)  -  0.10 
C 

DO  140  1=  M.N 

140  Till  rA(l)  =  STAR  I  +  DFi;n  Vl-M) 

C 

DO  160  K  =  M,N 

ir(K.FQ.1.0R.K.FQ.NP)GO  TO  llOR 
AI  I  =  (11  II:TA(K)  +  T1  IETA(K- 1  ))/2  0 
IJLI  =  (TI1F1  A(K)  +  T11ETA(K  +  1))  2  0 

C  GW  AND  GM  ARE  1  WO  ADJUS1 ABI  F.  PARAMFU  RS  DI  NOUNG  KIF  WID1 11  OF  TIIF 
C  GUASSIAN  PEAK  AND  THE  POWER  FOR  I  IIF  GUASSIAN  DIS I  Rini'TION  RFSPFCTIVFFI 
1108  GW  =  0.01 
GM  =  2.5 
C 

DO  150  J  =  M.N 
TTJ  =  (TT(J)  +  TT(J  +  l))/2.0 
D1FE  =  THETA(K)-TTJ 

YNORM  =  (1.0/(SQRT(GM*PI)*GW))*(GAMMA(GM)/GAMMA(GM-0.5)) 

DR  =  YNORM/(1.0  +  (DIFF**2)/(GM*(CiW**2)))**GM 
IF(K.FQ.  I.OR.K  EQ.  NP)GO  TO  150 

IF(1T(J),GF.ALI  AND.  TT(J)  ET.IM l)l’nP(K)  =  IIBP(K)  +  P{I) 

150  A1NTFN(K)  =  AINTEN(K)+  P(J)*BR 
160  CONTINLl 
C 

C  TMF  I  (HI.OWINC;  ARE  THE  WRITE  SI  AI  FMENTS  IN  THE  PROGRAM 
C 

WRITE(I6.I) 

DO  15  I  =  l.NX 
WRI11',(  16,.5)X2(I).SX(I).SY(I) 
l.s  CONTINF'I 
WR1TF.(26,2) 

DO  125  1=  l.NX 
EX(I)  =  FX(I)*  1000.0 
EY(1)  =  EY(I)*  1000.0 
EZ(I)  =  E7(1)*  1000.0 
125  CONTINUE 
DO  25  1=  1,NX 

WRriIX26,6)X2(l),FX(l).FY(l),F.7.(l) 

25  CONTINUF 
WRITE(36..3) 

DO  35  1=  1,NP 
x(n  =  xai  -  0.25 
WRITF(.5fi,7)D(I),X(I).IT(l) 
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35  CONTINIT 
WRHf;(46,4) 

DO  145  I  =  M,N 

WRrrE(46,8)  THF;rA(I).AINTF.N{I) 

145  CONTI  MJF 
WRITE(56,11) 

DO  55  1=  1.7 

WRITE(56,9)SPSI(r).D0S(I).D0R(I),D90S(I).D90F(I) 
55  CONTINUE 


STOP 

END 

C 

C 

C  THIS  ROUTINE  CAI  CULATES  1  Hi:  D-SPACING  AND  HIE  I  WO  rHi;i  A  VALUES  FROM  Ti 
C  STRAIN  EZ  1  HIS  ROUTINE  CAEDS  HIE  SUBROU  I INE  "SI  RAIN  l  O  GI  I  HIE  VAI  11 
C  OF  THE  STRAIN,  EZ.  THIS  ROUTINE  AESO  CAI.ES  ANO  I  HER  SUnROUHNl  'DP  I  S 
C  WHICH  IS  USED  TO  FIl  A  CURVE  THROUGH  HIF  POINTS  OF  D  VS  X  GENFRAH  D 
C  IN  SUBROU  I  INE  'DXIIT 

(^t^^t^^*****4****************»***********************»****************,* 

SUBROUTINE  DX2TT(D0.NX,EZ,D,  FI ) 

DIMENSION  EX(100).EY(100).F.Z(  100), D(2n00), 11(2000) 

C  IMPElCn  REAI  MIA-H.o-Z) 

PI  =  ACOS(-1.0) 

RAD  =  2.0H80.0/P1 
ALAMDA=  1.54178 
DO  =2.0235 
DO  100  1=  1,NX 
D(I)  =  (l.0  +  EZ(I))*D0 
1T(I)  =  ASIN(AEAMDA/(D(I)*2.0))*RAD 
100  CONHNUE 
RETURN 
END 
C 

Qt^^t******************************************»****r»****************t* 

C  THIS  ROU  FINE  CALCULATES  THE  D-SPACING  VALUES  AT  HIl-  NO  FCH  HP  FOR 
C  DIFFFRENT  VALUES  OF  PSI.  FOR  BOHl  SURFACE  AND  EDGE.  THF^SE  VALUES 
C  CAN  BE  GENERATED  FOR  ANY  POSIHON  ALONG  HIF  X-AXIS,  BUT  FOR 
C  VERIFICATION,  THE  VALUES  AT  HIF,  N01 CH  HP  ARF.  GIVEN  HERE 

(-*t****.*.*********************t****t****t***. ********. ************ . 

C 

SUBROUTINE  DVSSIN(inX2.EX,EY.EZ.SPSI,D0S,D0E,D90S.n90E.) 

DIMENSION  X2(100).EX(I0O),F.Y(IO0).F.Z(l00).SPSI(50) 

DIMENSION  D0S(50).D90S(5O).D0E(5O),D90E(50) 

PI  =  ACOS(-I.) 

PSI  =  0,0 
J=  I 

DO  20  1=  1,7 
PHI  =  0.0 

SPHI  =  SIN(PHI)«*2 
CPHI  =  COS(PHI)**2 

P0SH;  =  D0*(i,O-f-F.Z(J)  +  (EX(J)*CPHI  +  FY(J)*SPHI-EZ(.I)) 

♦  ♦(SIN(PSI)**2)) 

DOE(I)  =  D01 1 .0  +  EX(.I)  +  (F,Y(J)*C  PHI  +  E.Z(J)*SPI II-EXI.l)) 
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♦  ♦(SIN(PSI)*»2)) 

Pin  =  900*PI/180. 

SPm  =  SIN(PHF)**2 

cpm  =  cos(pni)»*2 

D90S(I)  =  D0*(  1 .0  +  EZ(J)  +  (EX(J)*CPI  H  +  E  Y(J)*SPm-r.Z(J)) 

♦  ♦(SIN(PS1)*»2)) 

D90E(I)  =  D0*(  1 .0  +  EX(J)  +  (F.Y(J)*C:PHI  +  EZ(J)*SPIII-FX(J)) 

♦  *(S1N(PS1)**2)) 

SPSI(I)  =  S1N(PSI)**2 
PSI=PS1+  1 5.0*  ACOS(- 1.0)/ 180  0 

20  CONTINUE 
RETURN 
END 
C 

Q^t^t^t***************************************************************** 

C  mis  SUBROUl  IN’E  CALCUEA7ES  Till:  STRESS  AND  S  ERAIN  VAEUES  IN'  THE  X 
C  AND  Y  DIRECTION  AS  A  EUNCTION  OE  THE  DIS  l  ANCi;  BE  I  WFT  N  NOTCHES 

C:  THE  S  ERAIN  VAEUES,  7  HE  TOTAI,  NE'MBI-R  ()»'  POIN  I  S  A1  ONG  EHi;  X  DIRI  C- 

C:  TION  AND  THE  CORRESPONDING  X  AND  Y  VAEUES  ARE  THEN  RETURNED  1  (1 

C:  THE  MAIN  PROGRAM  EOR  GENI  RAEING  THE  D  SPACING  PROEM  E 

Q't^4*********************t*»*****************t************************t* 

c 

SUBROE' I INE  STRAIN(Yl,X2.EX,EV.i:Z.SX.SY,NX) 

C  IMPEICIT  REAE*8(A-H.O-Z) 

DIMENSION  CX(100).EY(100),EZ(IOO), 

1  SX{  100),S  Y(  I00),TX  Y(100),X  I(  100),X2(  100) 

F.=  10300.0 
PNU=  1.0/3.0 
T  =  0.342 
A  =  0.250 
RHO  =  0.125 

V0=  A1  AN(SQRT(A/RHO)) 

C 

C  E  IS  THE  YOUNG  S  MODULUS  OE  HIE,  MATERIAL 
C  PNU  IS  THE  POISSON'S  RATIO. 

C  T  IS  THE  THICKNESS  OF  THE  SPECIMEN. 

C  A  IS  HALE  THE  DISTANCE  BETWEEN  NOTCHES 
C  RHO  IS  THE  NOTCH  ROOT  RADIUS  OF  CURVATURE 

c  VO  IS  mi:  hyperbola  parameter. 

C 

S  =  70.0 

c:  S=P/(2*A*I) 

C  SA  =  35.0 
Cl  =SIN(V0) 

C2=COS(V0, 

CA  =  S*CI/(V0+CI*C2) 

C 

C  S  IS  THE  NOMINAL  STRESS  EOR  THE,  LOAD  P. 

C  SA  IS  HIE  STRESS  AMPLITUDE  WHICH  IS  S/2. 

C  CA  iS  A  CONSS  I  ANT  WITH  THE.  DIMENSIONS  OE  STRESS 

C  C  IS  A  SCALING  I  AC'  EOR  USED  IN  THE  STRESS  DISTRIBl  HON  I  Ql  AHONS 
C  VI  =  SlNH(U)*COS(V) 

C  Xl  =  X/C  =  COSH(U)*SFN(V) 

C  U  AND  V  ARE  THE  CURVILINEAR  COORDINAI  ES 
C 
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C  =  A/C1 

C  THF  FOLl  OWING  IF  LOOP  CA1-CIJI,A'I  FS  TIIF  SI  RFSS  VAFUFS  IN  TUF  X  AND 
C  Y  DIRFCI  IONS,  SX  AND  SY.  SHEAR  S PRESS  TXY  AND  Till-  CORRF^SPONDING 
C  STRAIN  VALUES  EX,  EY  AND  EZ  AS  A  FUNCl  ION  OF  HIE  LINE  ALONG  THE 
C  BASE  OF  THE  NOTCHES.  THE  VAI.UES  ARE  FIRST  DERIVED  IN  CURVILINEAR 
C  COORDINATES  AND  A  COORDINATE  TRANSFORMAl  ION  IS  USED  TO  CONVERT  TO 
C  CARTESIAN  COORDINATES. 

C  X-MAX  IS  THE  MAXIMUM  VALUE  ON  THE  HYPERBOLA 
C  X-MIN  AND  XX-MIN  ARE  THE  MINIMUM  VALUES  ON  Fill:  HYPERBOI  A 
C  H  IS  THE  FACTOR  OF  DISTORTION  FOR  THE  HYPERBOLA 

C 

Ul=ALOG(Yl/C2  +  SQRT((Yl/C2)**2+  I)) 

XXMIN  =  -C*COSH(UI)*SIN(VO) 

XMAX  =  C*COSH(UI)*SIN(VO) 

XMIN  =  XXMIN 
NX  =  0 

X2(1)=XMIN 
J  =  0 

20  J  =  J  +  1 

XI(J)  =  XMIN/C 
Ql=  1.0  +  (XI(J)**2)  +  (YU*2) 

Q  =  SQRI  ((0  l-SQRT(OI  **2-4.0*(X  i(J)*‘2)))/2.n) 

V  =  ASIN(Q) 

U  =  ALOG(  Y  l/COS(V)  +  SQRT((Y  l/COS(V))»»2  +  1 .0)) 

II  =  S0RT((Y1  ♦*2)/COS(V)*»2  +  COS(V)**2) 

SU  =  CA*COSH(U)*COS(V)*(2.0*(H*»2)  +  COS(VO)**2-COS(V)**2)/II*M 
SV  =  CA*COSH(U)*COS(V)*(COS(V)**2-COS(VO)*»2)/II*M 
TUV  =  CA*SINH(U)*SIN(V)*(COS(VO)»*2-COS(V)**2)/IDM 
ALPHA  =  ACOS(C/IDCOSH(U)*COS(V)) 

SX(J)  =  (SU  +  SV)/2.0  +  (SU-SV)*COS(2*ALPI  I  A)/2.0-TU  V»SIN(2*  ALPI I  A) 

SY(J)  =  (SU  +  SV)/2.0-(SU-SV)*COS(2*ALPHA)/2.0  +  TUV*SIN{2*  ALPHA) 

TX  Y(J)  =  (SU-SV)*SIN(2*ALPHA)/2.0  +  TUV*COS(2‘ALPHA) 

THE  FOLLOWING  ARE  THE  CORRESPONDING  ELASTIC  S  PRAINS. 

EX(J)  =  (SX(J)-PNU*S  Y(J))/E 
EY(J)  =  (SY(J)-PNU*SX(J))/E 
EZ(J)-  -(SX(.I)  +  SY(J))*PNU/E 
C  S 1  =  (SX(J)  +  SY(J))/2,0  +  SQRT(((SX(.I)-SY(J))»*2)/4.n  +  PX  Y(J)*  *2) 

C  SI  =(SX(.I)+  SY(J))/2.0-SQRT(((SX(J)-SY(J))**2)/4n  +  TXY(J)**2) 

NX  =  NX  +  I 

XMIN  =  XMIN  +  (XMAX-XXMIN)/50.0 

K  =  J+  1 

X2(K)  =  XMIN 

IF(XMIN.I  F..0.26)GO  TO  20 

RETURN 

END 
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Sample  input  files 


C  Input  file  for  (200)  plane 
C 

2.0235 
0.00 
-0.5 
0.067 
2.01441 
-0.175 
0.0435 
2.01483 
-0.063 
0.0254 
2.01556 
-0.03333333 
0.01633.333 
2.01625 
-0.025 
0,0125 
2,0166875 
C 

C  Input  flic  for  (400)  plane 
C 

1.01175 
0.00 
-0.25 
0.0335 
1.00721 
-0.0625 
0.01875 
1.0075 
-0.03.3333.33 
0.013 
1.0077733 
-0.01666667 
0.00817 
1.0081229 
-0.01 
0.0053 
1 .00843 
C 

C  Input  flic  for  {.3.3.3)  plane 
C 

0.7788455 

0.00 

-0.1731 

0.025 

0.775.34771 

-0.0656 

0.0165 

0.7755154 

-0.0259 

0.01 


«) 

(y.) 

(ai) 

(bi) 

(Cl) 
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0.77578.-^1 

-0.0138 

0.006632 

0.7760245 

-0.011 

0.005495 

0.7761.397 
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Output  generated 


THF,  OUTPUT  THAT  CAN  BF.  GENFRATFD  FROM  TUI-  PROGRAM  ARF 

1)  STRESS  VALUES  IN  THE  X  &  Y  DIRECTIONS  VARYING  A1X)NG  X 
THESE  CAN  BE  OBTAINED  FOR  Dll  FEREN  I  HEIGHTS  ALONG  Y 

2)  STRAINS  IN  THF  X.  Y  AND  Z  DIRECTIONS. 

3)  D-SPACING  VALUES  CORRESPONDING  1 0  THF  STRAIN  FZ 

4)  BROADENED  INTENSITY  PROFILE  FROM  I  HE  D-SPACING  GRADIENT. 

5)  INTENSITY  PROFILES  FROM  TEN  SEGMENTS  DUE.  1  O  PARTIAL 
SUMMATION 
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